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SUMMARY 


The  state-of-information  is  quite  different  for  the  various  effects 
that  che  atmosphere  has  on  the  propagation  of  laser  radiation.  Molecular 
line  absorption  is  quite  well  understood,  and  several  sophisticated  line- 
by-line  computer  codes  and  comprehensive  line  parameter  compilations  are 
available.  On  che  other  liand,  the  mechanism  for  the  continuum  molecular 
absorptions  is  still  the  sublect  of  controversy,  and  there  is  consider¬ 
able  disagreement  between  various  Measurements  and  models.  Molecular 
.carter mg  theory  is  well  established  and  supported  by  measurements  and 
models. 

Theories  and  models  for  aerosol  absorption  and  scattering  are  also 
well  established,  but  aerosol  characteristics  of  the  atmosphere  (partic¬ 
ularly  the  size  distribution)  are  highly  variable  and  difficult  to  char¬ 
acterize.  As  a  result,  che  accuracy  of  most  predictions  of  aerosol  ef¬ 
fects  is  uncertain,  unless  detailed  measurements  of  aerosol  size  distri¬ 
butions  are  available  for  the  situation  of  Interest.  No  better  than 
order-of-magnitude  predictions  should  be  expected  from  correlations  of 
aerosol  extinction  in  different  spectral  bands  (e.g. ,  the  correlation  of 
infrared  transmission  with  the  visibility). 

The  effects  of  weak  to  moderate  levels  of  atmospheric  turbulence 
on  laser  radiation  transmission  is  well  understood.  In  this  "linear 
region"  of  turbulence  effects,  theories  and  models  are  available  for 
parameters  such  as  irradiance  variances  and  others  of  less  interest 
here  (o.g.,  beam  spread  and  wander,  and  polarization  and  coherence  ef¬ 
fects).  A  generally  accepted  theory  is  not  available  for  nonlinear  ef¬ 
fects  associated  with  higher  levels  of  turbulence  in  the  "saturated"  or 
"supersaturated"  region.  However,. there  are  a  number  of  satisfactory 
empirical  models  for  this  situation.  Once  again,  there  .is  considerable 
uncertainty  in  the  accuracy  of  available  models  for  predicting  the  state 
of  the  atmosphere,  this  time  with  regard  to  the  expected  turbulence 
levels;  also,  measuring  the  level  of  atmospheric  turbulence  in  the  field. 
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is  difficult.  The  situation  in  this  area,  however,  is  not  so  serious 
since  the  effects  of  interest  are  not  so  sensitive  to  the  level  of  tur¬ 
bulence,  and  the  maximum  turbulence  effects  are  somewhat  more  predict¬ 
able  than  the  aerosol  effects. 

Briefly,  for  molecular  effects,  we  are  recommending  the  use  of  the 
Air  Force  Geophysics  Laboratory's  (AFGL)  computer  code  LASER  to  generate 
simpler,  user-oriented  algorithms  for  rapid  prediction  of  molecular  ex¬ 
tinction  for  various  path  conditions.  Such  algorithms  have  been  developed 
for  a  limited  number  of  laser  lines.  Besides  treating  molecular  absorp¬ 
tion  (line  and  continuum)  and  scattering,  the  code  LASER  also  treats  aero¬ 
sol  extinction.  However,  to  be  conservative  in  safety  considerations,  we 
recommend  that  the  minimum  aerosol  effects  be  considered  (i.e.,  the  Clear 
Model  in  LASER).  This  recommendation  is  made  because  aerosol  attenuation 
is  so  poorly  predicted,  from  the  easily  obtainable  atmospheric  parameters 
such  as  visibility,  humidity,  and  wind  speed,  using  currently  available 
models.  For  predicting  the  turbulence  condition  of  the  atmosphere,  we 
recommend  using  Hufnagel's  1978  analytical  model  with  a  small  adjustment 
of  parameters  to  bring  it  into  better  agreement  with  measured  data.  For 
predicting  the  effects  of  atmospheric  turbulence  on  irradiance  statistics, 
we  recommend  using  the  classical  Rytov  expression  in  the  linear  region  and 
an  empirical  formula  developed  by  Johnson,  et  al.,  for  the  saturated  and 
supersaturated  regions.  It  appears  that  the  most  turbulence  can  do,  99 
percent  of  the  time,  is  to  increase  the  local  value  of  irradiance  over 
its  average  value  by  a  factor  of  five. 
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ATMOSPHERIC  EFFECTS  UPON  LOW-POWER  LASER  BEAM 
PROPAGATION 


INTRODUCTION 

The  main  objectives  of  this  study  were  (1)  to  establish  the 
.  state-of-.inf  ormation  regarding  the  transmission  characteristics  of  the 
atmosphere  for  low-power  laser*  radiation,  and  (2)  to  recommend  models 
for  predicting  those  characteristics.  The  models  are  to  be  employed 
for  predicting  safe  ranges  for  personnel  in  the  vicinity  of  low-power 
lasers  and  for  operational  planning  that  addresses  the  associated  safety 
considerations.  The  study  was  to  consider  the  absorption  and  scattering 
effects  of  natural  atmospheric  molecules  and  aerosols,  as  well  as  the 
effects  of  atmospheric  turbulence  over  arbitrary  slant  paths  in  the  at¬ 
mosphere.  Of  primary  concern  for  safety  considerations  are  the  predic¬ 
tion  of  the  total  average  molecular  and  aerosol  extinction  and  the  ir- 
«  radiance  statistics  associated  with  turbulence.  Models  were  also  re¬ 
quired  to  treat  various  meteorological  conditions,  including  rain,  snow, 
haze,  fog,  various  types  of  clouds,  and  atmospheric  turbulence  at  vari¬ 
ous  altitudes.  Spectral  regions  of  interest  include  the  ultraviolet, 
visible,  and  infrared  (i.e.,  approximately  from  0.2  to  10  hm) ,  with  the 
emphasis  on  tho:-'e  regions  at  the  frequencies  of  commonly  used  lasers. 

The  first  six  months  of  this  study  were  devoted  to  preparing  a 
bibliography  that  was  intended  to  \>v.  is  comprehensive,  and  user  oriented 
as  possible.  In  this  report  the  major  features  of  this  bibliography 
will  be  reviewed,  our  overall  conclusions  and  recommendations  will  be 
summarized,  and  available  information  will  be  evaluated.  The  details 
of  our  conclusions  and  recommendations  regarding  the  best  models  for 
the  present  purposes  will  also  be  given. 


*Lasers  for  which  nonlinear  effects,  such  as  thermal  blooming  and 
aerosol  modification,  can  be  neglected. 
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AVAILABLE  INFORMATION 


A  comprehensive  bibliography  of  information  relating  to  atmospher¬ 
ic  effects  on  the  transmission  of  low-power  laser  radiation  was  compiled 
during  the  first  six  months  of  this  study  [ll.  This  bibliography  con¬ 
sidered  information  on  theories,  measurements,  and  models  for  the  ef¬ 
fects  of  atmospheric  molecules,  aerosols,  and  turbulence  on  the  propaga¬ 
tion  of  laser  radiation  in  the  ultraviolet,  visible,  and  infrared  spec¬ 
tral  regions.  In  peripheral  areas  such  as  the  condition  and  composi¬ 
tion  of  the  atmosphere  and  the  optical  properties  of  aerosol  materials, 
only  the  principal  papers  were  included;. 

A  number  of  sources  were  utilized  to  compile  the  information  pre¬ 
sented  in  the  bibliography,  including  two  different  compucerized 
searches,  several  earlier  technical  reviews,  and  discussions  with  ex¬ 
perts  in  the  field.  The  most  useful  source,  in  terms  of  the  number  of 
reports  located,  was  the  Lockheed  computerized  bibliographic  search 
service  called  DIALOG.  This  system  indexes  papers  from  approximately 
4400  technical  journals,  1000  conference  or  symposium  proceedings,  and 
selected  books  and  U.S.  Government  reports.  Based  on  user-supplied  key 
words  and  phrases,  the  system  locates  pertinent  papers  and  prints  the 
reference  material  (e.g.,  title,  author,  report  or  journal  number),  as 
well  as  an  abstract  when  one  is  available.  This  search  and  one  by  the 
National  Technical  Information  Service  (NTIS)  included  many  papers  that 
were  not  nerf.  inert  to  present  problem  (e.g.,  papers  primarily  on 
laser  fusion,  flow  i nostics,  nonlinear  effects).  Thus,  considerable 
manual  review  of  the  titles  and  abstracts  was  required;  less  than  20 
per',  l  oi  the  references  listed  in  che  two  computer  searches  were 
deemed  pertinent. 

Over  1100  citations  to  pertinent  reports  and  papers  are  presented 
in  this  bibliography,  along  with  approximately  50  additional  references 
to  texts.  Over  80  percent  of  these  citations  are  annotated  with  a  few 
sentences  describing  the  main  points  of  the  paper.  This  bibliography 
is  arranged  in  ten  categor.es  as  shown  in  Table  1.  This  table  also  in¬ 
dicates  the  breakdown  as  to  where  the  information  appears  and  how  many 
in  each  category  are  annotated. 
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TABLE  1.  BIBLIOGRAPHY  SUMMARY  -  ATMOSPHERIC  EFFECTS  ON  LASER 
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A  thorough  cross-reference  to  59  subcategories  is  also  included 
in  this  bibliography.  This  cross-reference  specifies  (by  citation  num¬ 
ber)  each  citation  that  pertains  to  each  subcategory,  with  any  single 
citation  possibly  appearing  in  several  different  categories.  These  sub¬ 
categories  include  items  such  as  report  type  (measurement  report,  theo¬ 
retical  report,  etc.);  various  types  of  molecular  absorption  (line,  con¬ 
tinuum,  and  those  due  to  various  species  such  as  1^0  and  CC^);  various 
kinds  of  aerosols  vdust,  fog,  rain,  etc.)  and  their  characteristics 
(optical  properties  and  size  distributions);  various  turbulence  effects 
and  conditions;  and  specifically  addressed  lasers.  Another  feature  of 
this  bibliography  is  the  inclusion  of  synopses  of  telephone  conversations 
with  several  key  researchers  in  this  field  regarding  their  latest  work 
and  future  plans. 


EVALUATION  OF  AVAILABLE  INFORMATION 


In  order  to  establish  the  state-of -knowledge  in  the  field  of  atmo¬ 
spheric  effects  on  the  transmission  of  laser  radiation  and  to  recommend 
models  for  predicting  these  effects,  it  was  necessary  to  evaluate  the 
available  information  discussed  in  the  previous  section.  In  this  eval¬ 
uation,  we  concentrated  on  areas  of  primary  importance  (i.e.,  atmospher¬ 
ic  effects  that  would  have  the  most  potential  impact  on  eye-safety  con¬ 
siderations)  and  on  the  more  controversial  areas  where  the  theories,  ex¬ 
periments,  and/or  models  are  uncertain.  Thus,  we  examined  in  detail  the 
more  important  papers  in  the  bibliography  (and  some  others),  and  compared 
and  evaluated  their  results  and  conclusions. 

Information  on  ihe  three  main  atmospheric  effects  of  interest — 
those  due  to  molecules,  aerosols,  and  turbulence — are  evaluated  separate¬ 
ly  in  the  following  three  sections.  Due  to  the  importance  of  and  contro¬ 
versy  about  these  effects,  the  evaluation  in  each  section  concentrates 
on  the  molecular  continuum  absorption,  the  comparison  of  aerosol  extinc¬ 
tion  models  with  measurement,  and  the  expected  state  of  turbulence  in 
the  atmosphere.  Each  section  also  discusses  the  generally  accepted  the¬ 
ory  underlying  each  main  atmospheric  effect  and  the  more  pertinent  ex¬ 
periments  and  models  pertaining  to  that  effect. 

Molecular  Absorption  and  Scattering 

Molecules  produce  the  most  ubiquitous  effect  that  the  atmosphere 
has  on  the  transmission  of  laser  radiation  since  molecules  are  by  defi¬ 
nition  always  present  in  the  atmospheric  path;  whereas  aerosols  and  tur¬ 
bulence  need  not  be  present  in  significant  amounts.  The  number  density 
of  these  molecules  decreases  rapidly  with  altitude  above  sea  level  (gen¬ 
erally  exponentially),  so  attenuation  at  high  altitudes  is  much  less 
than  at  low  altitudes.  Several  molecules  in  the  atmosphere  (naturally 
occurring  and  others)  strongly  attenuate  radiation  in  the  ultraviolet, 
visible,  and  infrared  spectral  regions.  The  relative  proportion  of  many 


II 


of  these  molecules  (e.g.,  CO^,  CO,  CH^)  is  quite  uniform  both  spatially 
and  temporally  throughout  the  atmosphere.  The  relative  concentration  of 
others,  notably  H^O  and  0^,  varies  strongly  with  locale  and  time  and  does 
not  follow  a  general  exponential  distribution  like  the  other  species. 
Thus,  to  determine  the  eye-safe  range  to  a  laser  operating  in  the  atmo¬ 
sphere,  knowing  the  concentration  of  several  pertinent  molecules,  and 
their  statistical  variability,  is  also  important. 

Molecules  affect  the  transmission  by  two  different  mechanisms: 
absorption  and  scattering.  The  absorption  effects  are  commonly  further 
separated  into  those  due  to  individual  absorption  lines  causing  fine 
spectral  structure  and  those  due  to  "continuum  absorption"  varying 
smoothly  with  wavelength.  The  mechanism  producing  continuum  absorption 
is  not  well  understood,  but  is  thought  to  be  caused  by  either  the  far 
wings  of  a  large  number  of  individual  lines  or  by  the  existence  of  dimers 
or  polymers  of  certain  molecules,  particularly  water  vapor.  The  mecha¬ 
nisms  of  both  resonant  line  absorption  and  scattering  are  well  understood 
and  explained  by  classical  wave  mechanics.  This  section  evaluates  the 
various  theories,  measurements,  and  models  relating  to  atmospheric  mo¬ 
lecular  concentration,  molecular  line  absorption,  molecular  continuum 
absorption,  and  molecular  scattering. 

Atmospheric  Molecular  Concentration — The  attenuation  of  laser  ra¬ 
diation  by  a  particular  molecular  species  in  the  atmosphere  is  generally 
proportional  to  the  number  density  of  that  species  in  the  path  of  the 
beam,  although  self-broadening  effects  in  continuum  absorption  involve 
a  quadratic  dependency.  As  mentioned  above,  the  concentration  of  some 
of  the  more  attenuating  molecules  can  be  highly  variable  with  time  and/ 
or  locale;  thus,  it  is  very  important  to  have  information  on  the  concen¬ 
tration  of  certain  pertinent  molecules  and  their  statistical  variability. 
The  molecules  of  primary  concern  are:  1^0,  0^,  CC^,  CH^,  C^,  CO,  and 
^0;  all  except  the  first  two  are  generally  considered  uniformly  mixed 
throughout  the  atmosphere  (i.e.,  uniform  relative  concentration). 

The  measurement  of  these  and  other  molecular  constituents  of  the 
atmosphere  has  been  the  subject  of  numerous  studies  attempting  to  char¬ 
acterize  their  variation  with  altitude,  locale,  time,  ate.  Although  the 
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literature  survey  (1]  performed  in  the  first  phase  of  this  study  did  not 
emphasize  this  topic,  35  references  on  it  were  given.  Of  particular  in¬ 
terest  in  recent  years  have  been  measurements  of  0^  concentration.  Water 
vapor  and  carbon  dioxide  have  historically  been  the  subject  of  many  mea¬ 
surement  programs. 

Most  previous  calculations  and  analyses  of  atmospheric  attenuation 
of  radiation  have  been  based  on  standard  models  for  molecular  concentra¬ 
tions  that  were  derived  by  the  Air  Force  Geophysics  Laboratory  (AFGL) 
[2,3]  from  measurements  such  as  these.  Values  from  these  models  are  re¬ 
produced  in  Table  2.  A  more  recent  model  of  atmospheric  concentrations 
[4]  considers  more  molecules  and  gives  a  range  for  their  expected  vari¬ 
ation  with  altitude,  latitude,  and  time  of  year  (see  Table  3).  It  is 
noted  that  there  are  large  variabilities  associated  with  the  concentra¬ 
tion  of  several  of  these  molecules.  Water  vapor  concentration  is  of 
particular  concern  because  it  is  so  highly  variable  (particularly  with 
altitude  and  latitude)  and  because  this  molecule  is  so  strongly  absorb¬ 
ing  in  several  spectral  regions.  Figure  1  shows  that  the  measured  vari¬ 
ation  in  water  vapor  concentration  witli  altitude  for  several  seasons 
does  not  correlate  well  with  the  AFGL  models.  These  models  generally 
indicate  a  much  slower  reduction  in  concentration  with  altitude,  which 
is  not  conservative  with  regard  to  eye  safety.  For  this  reason,  it  is 
advisable  to  base  attenuation  predictions  on  the  locally  (and  currently) 
measured  water  vapor  concentations  whenever  possible,  rather  than  on 
values  from  a  standard  model.  This  is  also  advisable  for  other  species, 
but  is  not  as  critical. 

Molecular  Line  Absorption — Frequently  the  dominant  attenuation  of 
laser  radiation  in  the  atmosphere  is  due  to  absorption  by  individual 
vibration-rotation  lines  of  atmospheric  molecules.  This  mechanism  is 
particularly  important  for  radiation  produced  by  gas  lasers  using  a 
molecule  present  in  the  atmosphere,  such  as  CC^.  In  these  cases,  the 
emission  and  absorption  are  caused  by  the  same  transition  (i.e.,  reso¬ 
nant  absorption),  so  the  laser  frequency  is  right  at  the  center  of  the 
absorption  line  where  the  absorption  is  greatest. 
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TABLE  2.  CONCENTRATIONS  OF  IMPORTANT  ABSORBING 
ATMOSPHERIC  MOLECULES  [2,3] 


Uniformly  mixed 


Nonun if ormly  mixed 


Sea  level  density 


Concentration 
(parts  per 
million3) 

Model 

atmosphere 

(g/m3) 

Molecules 

h2° 

°3 

co2 

330 

Tropical 

19 

5.6xl0"5 

n2o 

0.28 

Midlatitude 

14 

6 .OxlO-3 

CO 

0.075 

summer 

Midlatitude 

3.5 

6.0xl0-5 

ch4 

1.6 

winter 

Subarctic 

9.1 

4.9xl0-5 

°2 

2.10xl05 

summer 

Subarctic 

1.2 

4.1xl0-5 

aBy  volume 

in  dry  air. 

winter 

TABLE  3.  CONCENTRATION  OF  ABSORBING  ATMOSPHERIC 
MOLECULES  AND  THEIR  VARIABILITY  [4] 


Average 


Molecule 

Average  concentration 
(volume  basis) 

Molecule 

concentration 
(volume  basis) 

N2 

0.78084 

NO 

io~8-.io~6 

°2 

0.20946 

no2 

io"9-io"6 

h2o 

1.3xl0-7  to  4 . 5xlO-2 

hno3 

2.8xl0~9 

co2 

3.10X10-4 

nh3 

<10~6 

CH  4 

(1-I.4)xl0-6 

so2 

(0.5-7.2)xl0-9 

h2 

5x1 0~7 

h2s 

(1.6-16)xl0~9 

CO 

(0.5-2.5)xl0~7 

hcho 

<io'7 

°3 

(2-7)xlO-8 

HC1 

(l-2.6)xl0~9 

N20 

(2.7-3.5)xl0“7 

N03,0H,H02,CH30 

5xl0-11 
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AFGL  MODELS 


Figure  1.  Atmospheric  water  content  aa  a  function  of  altitude  [5]. 


For  atmospheric  conditions  at  low  altitudes  (i.e.,  near  standard 
pressure),  the  laser  absorption  by  an  individual  absorption  line  is 
generally  well  expressed  by  the  simple  Lorentz  expression: 


k  =* 


nSy/ir 


(1) 


where 


k  =  Absorption  coefficient  per  unit  length;  i.e.  transmission 
T  =  exp(-kL),  where  L  is  the  path  length 

=  Laser  frequency 

n  ■*  Molecular  concentration  of  the  species  being  considered 
S  ■  Absorption  line  parameter  of  strength 
Y  *  Absorption  line  parameter  of  halfwidth 
Vq  **  Absorption  line  parameter  of  frequency. 


The  strength  depends  on  temperature  (T)  through  the  vibration  and  rota¬ 
tion  partition  functions  (Qy  and  Q^)  and  the  Boltzmann  factor  [3]: 


S(T) 


Mid  Mi)' 

’l.439(T-Ts  )E" 

QyCT)  Qr(T) 

exp 

IT 

s 

(2) 


where  Tg  is  some  standard  temperature  when  S  is  known,  and  E"  is  the 
energy  of  the  lower  state  of  the  transition.  The  halfwidth  is  propor¬ 
tional  to  the  local  atmospheric  pressure  (p),  and  generally  Inversely 
proportional  to  temperature  to  some  power: 


Y(T,P)  •  Y(r,.p,)(p/p.)(t,/l)"  (3) 

where  pg  is  a  standard  pressure  at  a  condition  where  Y  is  known.  The 
power  n  is  usually  taken  to  be  1/2  --  corresponding  to  temperature- 
independent  collision  diameters,  although  there  is  some  uncertainty  in 
its  actual  value. 
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The  broadening  of  an'y  particular  transitipn  results  from  inter¬ 
actions  with  all  of  the  molecules  of  the  gas  mixture.  The  separation 
of  the  contributions  from  the  different  molecules  is  a  tedious  task., 
though  fortunately  the  halfwidt.is  included  in  the  line  tabulation  make 
thij  separation  largely  unnecessary.  The  halfwidths  in  the  tabulation 
are  diluted  air  halfwidths,  characteristic  of  the  line  as  it  appears  in 
air  of  normal  composition.  However,  because  the  halfwidth  is  made  up  of 
so  many  contributions,  and  because  the  temperature  exponents  of  these 
contributions  vary  widely,  the  use  of  a  single  pressure  exponent  of  one- 
half  is  very  much  an  oversimplification;  so  large  temperature  correc¬ 
tions  should  be  avoided. 

The  AFGL  Line  Parameter  Compilation  (3]  provides  all  of  these 
parameters  for  the  most  important  atmospheric  absorption  lines.  To 
evaluate  the  total  absorption  due  to  all  lines  that  significantly  af¬ 
fect  the  laser  radiation,  it  is  necessary  to  sun  the  individual  absorp¬ 
tion  coefficients  from  all  of  these  lines.  Generally,  this  is  done  by 
considering  all  of  the  important  absorption  lines  within  a  certain  cut¬ 
off  frequency  of  the  laser  frequency  (usually  20-50  cm  ^) .  The  impor¬ 
tance  of  each  absorption  line  is  determined  by  the  absorption  It  would 
produce  at  its  center  for  an  extreme  atmospheric  path  tangent  to  the 
Earth's  surface  and  extending  from  space  to  space;  lines  producing  less 
than  10  percent  absorption  for  this  path  are  neglected  in  the  AFCL  com¬ 
pilation.  However,  this  cutoff  was  not  employed  in  two  situations: 

(1)  in  regions  of  very  strong  absorption  where  relatively  weak  lines 
above  this  limit  are  ignored;  and  (2)  for  Q-branch  lines  below  this 
limit  which  are  included  when  it  is  felt  that  their  cumulative  effect 
might  be  significant. 

At  higher  altitudes,  where  the  pressure  Is  too  low  for  collisions 
to  have  significant  broadening  effects,  the  Doppler  broadening  becomes 
dominant  and  the  absorption  coefficient  is  better  approximated  by: 


where  the  Doppler  halfwidth  is  expressed  as 


YD  -  —  [(2  £n  2)2kT/m],S 


(5) 


and 


y  -  (In  2)**  (\>L  -  v0)/yd  ( 

where  c,  k,  and  m  are  the  speed  of  light,  Boltzmann's  constant,  and 
molecular  mass,  respectively. 

For  intermediate  altitudes  (i.e.,  pressures),  where  the  Doppler 
and  pressure  broadening  are  of  the  same  order,  a  good  approximation  is 
obtained  by  convolving  the  two  profiles  to  get  what  is  referred  to  as 
the  Voigt  profile: 


k 


(yko/w) 


_ e*p(-t2) 

y‘  +  (x  -  t)2 


dt 


where 


k 

y 


o 


S/V 


Y/y, 


1 

2 


x  = 

yd 


(7) 


Other  more  complicated  line  shape  models  have  been  developed;  for  example 
[6],  the  full  Lorentz: 
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?  <VV,> 


_jc _ _ 

V#)2  +  Yr 


(vL  +  vj)a 


(8) 
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and  Van  Vleck-Weisskopf : 


i  <VV°>! 


_ y  t  Y  1 

-  v0)2  -f  Y2  +  +  Vo)2  +  Y2  I  ; 


(9) 


and  the  kinetic  model*: 


k  = 


4Sv  v0 

L 

7T 


[(yl  -  +  4vVj 


(10) 


but  generally  the  simple  Lorentz  and  Voigt  models  give  satisfactory 
agreement  with  measurements. 


Mainly  two  different  kinds  of  measurements  have  been  made  to  deter¬ 
mine  the  molecular  line  absorption  of  laser  radiation."**  The  first, 
direct  measurement  of  the  attenuation  of  a  laser  beam  through  a  path  of 
known  composition,  is  easiest,  but  does  not  give  enough  information  to 
extrapolate  to  other  conditions  with  confidence.  The  second  determines 
the  basic  line  parameters  (i.e.,  strength,  width,  and  position),  thus 
allowing  accurate  prediction  of  the  molecular  line  absorption  in  other 
situations.  Measurements  of  the  second  type  of  measurement  by  Woods 
et  al.  [8,9]  have  indicated  several  inaccuracies  in  the  original  AFGL 
Line  Parameter  Compilation  (see  Figure  2).  These  measurements  were  made 
on  air-broadened  samples  of  N20,  H DO ,  and  CH^  in  the  DF  laser  region  [8] 
and  on  air-broadened  H^O  samples  in  the  CO  laser  region  [9].  They  in¬ 
dicated  that 

(1)  agreement  was  excellent  for  the  N^O  line  parameters 

(2)  HDO  lines  tended  to  be  stronger  and  wider  than  calculated 
(by  factors  of  50  percent  and  25  percent,  respectively) 

(3)  measured  ClI^  spectra  were  many  times  stronger  than 
calculated  and  contained  many  additional  lines 

(4)  Loren tz  shape  may  be  more  nearly  correct  than  previously 
thought,  if  accurate  positions,  widths,  and  strengths 
are  used  and  If  the  wing  cutoff  is  not  restricted  to  the 
usual  25  cm-1. 

*Somu times  this  model  is  attributed  to  Zhevakin  and  Naumov  [7], 

**The  authors  wish  to  acknowledge  the  helpful  discussion  of  this 
subject  w I  Lb  I’.  Smith  of  OptlMetrJ.es,  Ann  Arbor,  Michigan. 
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TRANSMISSION  (%) 


Figure  2.  Selected  example  of  measured  and  calculated  spectra 
in  the  vicinity  of  the  P2(6)  DF  laser  line  [8]. 


Measurements  of  this  type  were  also  performed  in  the  HF  spectral  region 
iiu]  and  indicated  similar  results.  The  conclusions  of  these  latter 
two  references  i 9 , 1 0  j  were:  (1)  it  is  necessary  to  consider  absorption 
lines  in  this  region  as  far  as  200  cm  ^  from  the  laser;  and  (2)  the 
exponent  on  the  term  in  the  denominator  of  the  l.orcnt z  formula 

should  be  approximately  1.9  Instead  of  the  usual  value  of  2.  The  re¬ 
sults  of  the  latter  measurements  seemed  to  agree  better  (especially  for 
the  weaker  absorption  lines)  with  the  line  strengths  given  by  Flaud  and 
Camy-Peyret  [11]  than  with  those  of  Benedict  which  are  used  In  the  AFC.L 
compilation.  The  Naval  Kesearch  Laboratory  (NKL)  is  planning  similar 
measurements  using  a  tunable  laser,  but  the  results  will  not  be  avail¬ 
able  for  some  time. 
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There  have  been  survey-type  spectral  scans  [12]  in  spectral  re¬ 
gions  where  the  HBr  and  Xe  lasers  operate  that  indicate  no  apparent 
discrepancies  with  the  AFGL  Line  Parameter  Compilation;  however,  the 
survey  did  indicate  that  the  commonly  quoted  value  for  the  position  of 
the  Xe  laser  transition,  2145  cm  \  should  actually  be  2168  cm  .  This, 
of  course,  could  drastically  change  the  level  of  molecular  absorption 
calculated  by  models.  Many  discrepancies  have  been  accounted  for  in  sub¬ 
sequent  versions  of  the  compilation;  however,  only  a  small  part  of  the 
spectrum  (i.e.,  near  the  DF,  CO,  HF,  HBr,  and  Xe  laser  lines)  was  inves¬ 
tigated.  This  implies  that  problems  could  also  exist  in  other  spectral 
regions  where  lasers  operate. 

There  have  been  a  number  of  direct  laboratory  measurements  of  the 

attenuation  of  radiation  from  a  number  of  different  lasers  (i.e.,  CO2, 

CO,  DF,  HF,  and  erbium) .  These  measurements  were  made  primarily  in 

multiple-pass  White  cells  or  more  recently  with  spectrophones  under  a 

variety  of  experimental  conditions  (see  Table  4).  The  CO^  laser  is  the 

tnost  popular  for  this  kind  of  measurement.  McCoy  et  al .  [13]  made  one 

of  the  earliest  measurements  using  a  980-m  White  cell  with  a  mixture  of 

CO2  and  air.  They  scaled  their  measured  results  for  the  P20  line  to 

standard  atmospheric  conditions  (i.e.,  330  ppm  C0_  at  1  atm  total  pres- 

^  -1 

sure)  to  get  an  absorption  coefficient  of  0.0694  km  ,  compared  with 
their  calculated  value  of  0.076  km  ^  based  on  independent  measurements 
of  strength  and  width.  Henry  [18]  measured  an  intermediate  value  of 
0.073  km  The  discrepancy  in  these  basic  "dry  air"  measurements  and 
calculations  is  as  yet  unexplained.  Actually,  the  atmospheric  extinc¬ 
tion  of  CO2  laser  radiation  is  usually  dominated  by  the  water  vapor  con¬ 
tinuum  in  this  region,  as  will  be  discussed  in  the  following  section. 

A  number  of  field  measurements  of  the  extinction  of  CO2  laser  radiation 
have  also  been  made. 

Absorption  of  the  P20  CC^  laser  line  radiation  by  dry  air  mixtures 
has  been  measured  by  a  number  of  other  experimenters  (see  Table  5). 

There  seems  to  be  good  agreement  at  300  K  for  a  value  near  0.071  km"'1', 
and  Indications  of  a  reduction  in  absorption  with  reduction  in  tempera¬ 
ture.  Moskal  enko  et  nl  •  1 2 1  |  nJ.no  measured  the  absorption  of  the  P20 
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TABLE  4.  MAJOR  LABORATORY  MEASUREMENTS  OF  MOLECULAR  ABSORPTION  OF  LASER  RADIATION 


TABLE  5. 


l.AHOkATOKY  MEASUREMENTS  OF  DRY  AIR  ABSORPTION 
COEFFICIENT  FOR  THE  P20  CO  LASER 

Absorption 

Temperature  cool  f icient 


Exper imenter 

00 

Stephenson  et  al. 

29^ 

OppenheLm  and  Devir 

300 

Moskalenko  et  a). 

290-300 

Gerry  and  Leonard 

2/3 

McCu^bin  and  Mooney 

300 

McCubbin  et  al. 

300 

Long  and  McCoy 

- 

(kin'S  Reference 


0.075 

22 

0.071 

23 

0.055-0.071 

21 

0.040 

24 

0.071 

25 

0.073 

26 

0.029 

27 

CO.,  laser  line  in  water  vapor  and  showed  an  increase  in  absorption  with 
temperature  (i.c.,  absorption  coot t ic  icnts  of  0.055-0.071  km  ^  for  T  ” 
290-300  K,  respectively),  but  no  details  of  the  experiment  were  provided. 

Extensive  laboratory  measurements  have  also  been  made  of  the  molec¬ 
ular  absorption  of  OF  laser  radiation  by  atmospheric  gases.  Spencer 
et  al.  (In,  determined  the  ait  orpt imt  by  individual  samples  of  atmospher¬ 
ic  gases,  considering  gases  whose  line  attenuation  of  DF  radiation  is 
the  strongest  —  111)0,  N^O,  CH^ ,  and  CO^.  In  Figure  3,  the  total  attenu¬ 
ation  indicated  by  the  measurements  for  a  midlatitude  winter  model  is 
compared  wii.h  theoretical  calculations  using  the  AFCL  line  parameters. 

For  the  17  DF  lines  measured,  the  agreement  is  within  approximately  a 
factor  of  two.  It  should  he  noted,  however,  that  this  comparison  is  for 
one  of  the  drier  model  atmospheres  and  the  water  vapor  continuum  is  not 
considered,  so  much  larger  uncerta  ml  in  the  molecular  attenuation 
will  be  encountered  in  practice. 

More  recently,  Mills  i20|  measured  t  lie  attenuation  of  eight  DF 
laser  lines  as  a  function  of  the  concentration  of  these  same  gases.  All 
combinations  of  laser  lines  and  gases  were  not  considered  since  some  of 
the  lines  are  not  significantly  affected  by  some  of  the  gases.  The 


Figure  3.  Comparison  of  experimental  and  theoretical  absorption 
coefficient  values  for  17  DF  laser  lines  for  the 
midlatitude  winter  model  atmosphere  (H_0  continuum 
not  included).  [23] 

attenuation  of  N^O  of  five  laser  lines  [i.e.  P^_2(6),  P3_2^)>  P3-2^’ 

?2  ,(10),  ?2  ,(11)]  was  measured  to  be  linear  over  concentrations  rang¬ 
ing  from  1-12  ppm  up  to  50-220  ppm.  Attenuation  coefficients  thus  in¬ 
dicated  are  compared  in  Table  6  with  calculations  based  on  the  AFGL  line 
parameters,  and  with  measurements  by  Spencer  et  al.  i  1 6  J  and  Deaton 
et  a] .  [19].  Mills'  measurements  agree  best  with  theory  (i.e.  generally 
about  1  percent,  but  -8.6  percent  for  one  line).  Although  these  measure¬ 
ments  were  made  at  higher  than  normal  atmospheric  concentrations,  the 
linear  behavior  indicates  that  self -broadening  effects  still  are  not  im¬ 
portant  and  that  the  results  can  be  scaled  accurately  to  lower  concentra¬ 
tions.  Deaton's  measurements  were  made  at  up  to  1000  times  normal  atmo¬ 
spheric  concentration  and  are  calibrated  to  Mills'  data.  Also  shown  in 

Table  6  are  similar  measurements  for  CH. .  The  differences  with  theory 

4 

are  much  larger  for  this  molecule,  but  the  measurements  are  more  self- 
consistent,  indicating  a  potential  difficulty  in  the  values  contained  in 
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TABLE  6.  MEASURED  AND  THEORETICAL  N2O  ABSORPTION  COEFFICIENTS 
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early  versions  of  the  AFGL  line  compilation.  However,  the  absorption 

-2  -1 

coefficients  of  Table  6  are  typically  less  than  10  km  ;  hence  these 
coefficients  are  very  small  compared  with  those  due  to  other  atmospheric 
species. 

Mills  also  measured  the  attenuation  of  the  ^(8)  DF  laser  line 

radiation  due  to  pure  CO  at  pressures  of  248,  503,  and  761  torr  over  a 

^  -4-1 

0.7317-km  path  and  found  a  value  of  5.3x10  km  ,  assuming  a  self¬ 
broadening  coefficient  of  unity.  He  also  presented  unpublished  measure¬ 
ments  by  Meyers  of  General  Dynamic s/Convair  (see  Table  7)  that  indicated 
absorption  coefficients  of  some  lines  that  are  four  orders  of  magnitude 
higher  than  theory,  apparently  due  to  the  omission  of  an  isotopic  or 
weak  CO^  band  from  the  AFGL  compilation.  In  any  event,  this  absorption 
is  still  small  compared  with  that  due  to  H2O  continuum  or  HDO  lines. 

Mills'  measurements  of  HDO  absorption  of  six  D7  laser  lines  were 
fit  with  a  model  of  the  form 


k  =  ap  +  bp^  (11) 

where  k  is  the  absorption  coefficient  and  p  is  the  partial  pressure  of 
HDO-enriched  water  vapor.  These  coefficients  were  then  adjusted  for 
normal  isotropic  abundances  of  HDO  (i.e.,  0.03  percent)  and  the  expres¬ 
sions  evaluated  for  the  Midlatitude  Summer  Atmospheric  Model  (see  Table 
8).  The  differences  with  theory  (i.e.,  calculations  with  AFGL  line 
parameters)  are  seen  to  vary  between  30  and  142  percent.  These  measure¬ 
ments  are  subject  to  errors  due  to  unaccounted  effects  of  D20  as  well  as 
possible  absorption  or  condensation  of  water  on  the  mirrors  and/or  win¬ 
dows  of  the  White  cell. 

Measurements  of  the  absorption  of  CO  laser  radiation  by  atmospher¬ 
ic  gases  are  limited.  The  most  extensive  measurements  have  been  by  Long 
[17;  28-33]  using  ^0  vapor  broadened  by  ^  in  a  12-m  multipass  White 
cell.  Long  also  made  measurements  on  two  of  the  more  highly  absorbed 
lines  [32].  His  earliest  measurements  at  a  HjO  partial  pressure  of  8.89 
torr  [30]  had  considerable  scatter,  and  when  compared  with  theory  (i.e., 
using  'FGL  line  parameters),  the  mean  appeared  to  vary  from  a  factor  of 
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TABLE  7.  MEASURED  C02  ABSORPTION  COEFFICIENTS  FOR  SIX 
DF  LASER  LINES  (TAKEN  FROM  MILLS  [20]) 

.  ,  Absorption  coefficient  (km  *) 

Laser  frequency  - <- - — — 

DF  line  (cm-1)  Meyers  Mills 


P2-i(6) 

2680.179 

0.19  x 

10~% 

- 

P2-i(7) 

2655.863 

X 

00 

• 

10~* 

- 

P2-i (8) 

2631.068 

9.26  x 

lO'" 

5.3  x  nr* 

Ps-2(6) 

2594.198 

6.85  x 

io-** 

- 

Pj-2(7) 

2570.522 

5.24  x 

io-1* 

- 

00 

w 

w 

1 

As 

2546.375 

4.21  x 

10"% 

— 

TABLE  8.  HDO  ABSORPTION  COEFFICIENTS  EXTRAPOLATED 
TO  0.03  PERCENT  RELATIVE  HDO  ABUNDANCE 

Absorption  coefficient  (km-1) 
Measured 


Theory 


Line 

0.037.  HDO 

14. 

26  torr  H20 

14.26  torr 

P2-l (6) 

3.39 

X 

10~  3  p 

+ 

2.05 

X 

io~V 

5.24 

X 

10-* 

(38)a 

3.79 

X 

10~2 

P2-i (7) 

5.94 

X 

10~3p 

+ 

5.24 

X 

io-y 

9.54 

X 

io-1 

(30) 

7.35 

X 

io-* 

P2-i<8) 

3.38 

X 

10-l*p 

+ 

5.78 

X 

io-v 

6.00 

X 

io-3 

(-34) 

9.12 

X 

io- * 

P 3  —  2  (6) 

1.11 

X 

10-3p 

+ 

7.62 

X 

io-y 

1.74 

X 

io-1 

(142) 

7.18 

X 

10*  3 

Pj-2 (7) 

4.71 

X 

10-"p 

+ 

9.31 

X 

io-v 

8.61 

X 

10"  * 

(90) 

4.53 

X 

io- 3 

P  3  —  2 (8) 

1.45 

X 

10_,,P 

+ 

1.93 

X 

io-‘P* 

2.46 

X 

10"  3 

(118) 

1.13 

X 

10"  3 

(  )  Indicates  percentage  difference  from  theory. 
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1.15  high  at  total  pressures  of  126  torr  to  a  factor  of  0.65  low  at  767 
torr  (see  Figure  4).  His  next  measurements  (31]  were  not  compared  with 
theory,  but  he  did  state  that  the  absorption  was  considerably  higher  than 
predicted,  in  agreement  with  spectral  measurements  by  Woods  et  al.  [9] 
that  indicate  "super-Lorentz"  behavior.  The  reason  Long's  earliest  mea¬ 
surements  [30]  did  not  exhibit  this  behavior  (see  Figure  4)  is  thought 
to  be  tied  to  his  selection  of  laser  lines.  That  is,  all  of  the  lines 
were  very  close  to  the  center  of  strong  absorption  lines  (v  s  vT),  so 
the  different  exponent  on  the  (v-vq)  term  is  not  significant  for  these 
lines.  The  next  measurements  Long  reported  [31]  had  much  less  scatter, 
but  also  displayed  this  trend.  These  were  also  compared  to  a  super- 
Lorentz  shape  and  found  generally  in  good  agreement.  Long's  latest 
measurements  [32]  were  on  two  highly  absorbed  lines:  P^g(lO)  and  P^(12). 
As  expected,  the  latter  agreed  well  with  the  standard  Lorentz  model  since 
the  laser  line  is  near  a  strong  absorption  line  (see  Figure  5),  while  the 
former  required  the  super-Lorentz  model  because  it  is  farther  from  any 
strong  lines.  Also  shown  in  Figure  5  are  spectrophone  measurements  by 
Long  [32]  and  measurements  by  Rice  [34]  using  a  short  path  cell  (71  cm). 
Rice's  measured  absorption  coefficients  are  below  those  of  Long  and 
Lorentz  theory,  and  have  considerable  scatter.  Rice's  earlier  measure¬ 
ments  [35]  have  even  more  scatter,  but  generally  lie  above  standard 
Lorentz  theory  (see  Figure  6). 

A  number  of  experimenters  have  made  measurements  of  the  attenua¬ 
tion  of  laser-radiation  in  the  natural  atmosphere.  Generally  these 
measurements  are  not  as  useful  as  laboratory  measurements  because  of 
uncertainties  in  the  concentrations  of  the  molecules  and  aerosols  along 
the  path.  However,  field  measurements  can,  in  certain  cases,  add  to  the 
confidence  of  the  theory  and  laboratory  measurements,  or  possibly  indi¬ 
cate  shortcomings  in  them. 

An  example  of  field  measurements  lending  credence  to  laboratory 
data  is  shown  in  Figure  7,  The  solid  curve  representing  McCoy's  labo¬ 
ratory  measurements  [13]  of  CO^  laser  absorption  (at  330  ppm  CO^)  are 
slightly  below  most  of  the  field  measurements,  so  consideration  of  aero¬ 
sol  attenuation  could  bring  them  into  agreement.  The  field  data  from 


28 


H,0  PARTIAL  PRESSURE  -  8.89  ton 


absorption  coefficients  by 
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Figure  5.  Measured  and  calculated  data  for  water  vapor 
absorption  coefficient  at  1837*436  ci  132} 
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6.  Comparison  of  measured  and  Lorentz* calculated 
absorption  coefficients  by  Rice  (35]. 


WATER  VAPOR  PRESSURE  (torr) 


Figure  7.  Comparison  of  laboratory  and  outdoor  extinction 
data  for  P20  CO2  laser  line  (10.59  urn)  [4], 

McCoy  et  al.  [131  and  unpublished  data  by  Gilmartin  of  MIT  Lincoln 
Laboratory  are  in  very  good  agreement  with  the  laboratory  measurements 
and  tend  to  verify  the  nonlinear  dependence  on  water  vapor  content. 
Unfortunately,  the  source  of  Figure  7  [4]  did  not  give  complete  refer¬ 
ences  for  the  data  of  Goodwin  or  of  the  Naval  Research  Laboratory;  the 
data  attributed  to  Rensch  is  apparently  a  partial  set  of  that  given  by 
McCoy  et  al.  [  1  3  j . 

Extensive  outdoor  measurements  of  the  attenuation  of  DF  laser  ra¬ 
diation  were  conducted  by  NRL  [36  [along  a  5-km  path  at  the  Capistrano 
test  site  in  California  during  June  to  September  1975.  In  Table  9, 
these  measurements  are  compared  with  values  calculated  by  using  AFGL 
line  parameters  for  22  DF  laser  lines.  It  is  seen  that  the  measurements 
are  usually  about  the  same  as  the  calculations,  or  larger  (allowing 
"room"  for  aerosol  attenuation),  except  for  a  few  lines  [i.e.,  ?2(8), 
P1(9),  p2(4),  P1(7),  P^(6),  and  P^(5)].  Dowling  [36]  does  not  discuss 
this  problem,  but  it  is  interesting  to  note  that  these  lines  all  are 
near  the  larger  frequency  half  of  the  spectral  range  of  the  measurements. 

Later  experiments  by  NRL  [37],  at  Cape  Canaveral  Air  Force  Station 
during  the  spring  of  1977,  measured  the  atmospheric  attenuation  of  ra¬ 
diation  from  UeHe,  ild-YAG,  DF,  CO,  and  CO2  lasers.  These  measurements. 
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TABLE  9.  COMPARISON  OF  RADC  FIELD  MEASUREMENTS  WITH  LABORATORY  MEASUREMENTS 
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over  a  5.1-km  path  near  the  ocean,  were  thought  to  be  influenced  much 
n»re  strongly  by  aerosols  than  were  the  California  data,  probably  due 
to  changes  in  offshore  and  ocean  wind  conditions  at  the  Florida  site. 
Analysis  of  this  data  is  not  yet  complete. 

Another  well-controlled  measurement  of  the  atmospheric  attenuation 
of  the  radiation  of  four  DF  laser  lines  was  conducted  at  a  610-m-long 
outdoor  site  at  Rome  Air  Development  Center  [7].  This  site  is  approx¬ 
imately  24  km  from  the  nearest  urban  environment,  so  urban  and  industri¬ 
al  pollution  should  be  small.  However,  the  path  extends  0. 9-4.6  m  above 
low-lying  wet  grasslands  and  wet  swampy  areas,  so  the  water  content  in 
the  path  was  thought  to  be  higher  than  that  monitored  at  the  receiver 
and  transmitter  locations;  also,  high  levels  of  water-type  aerosols  were 
suggested.  These  field  measurements  are  compared  to  laboratory  measure¬ 
ments,  showing  that  they  are  factors  of  from  2  to  30  higher  than  labo¬ 
ratory  measurements.  In  Table  9,  this  factor  seems  to  correlate  well 
with  relative  humidity,  supporting  the  above  rationalization  of  these 
discrepancies.  The  laboratory  measurements  used  ir.  the  comparison  were 
a  composite  of  individual  measurements  on  various  atmospheric  species 
from  different  sources  —  Mills  {20],  Deaton  et  al.  [19],  Meyers  [ 38 J , 
and  Burch  et  al.  [39]  —  but  do  not  include  the  effects  of  aerosols. 
However,  even  after  accounting  for  nominal  aerosol  attenuation,  the 
field  measurements  are  still  generally  much  larger. 

Another  field  measurement,  by  Borisov  [40]  as  reported  by  Adiks 
et  al.  [41],  suffers  from  the  opposite  deficiency.  That  is,  measured  at¬ 
mospheric  attenuation  of  CO^  laser  radiation  seems  to  be  well  below  sev¬ 
eral  laboratory  measurements  f nd  theoretical  calculations.  This  implies 
difficulty  with  this  particular  field  measurement  by  Borisov,  since  the 
ether  data  are  in  general  agreement.  In  general,  however,  the  existing 
laboratory  and  field  measurements  tend  to  support  the  currently  accepted 
theoretical  information  on  molecular  lint*  absorption. 

Modeling  of  molecular  line  absorption  is  a  highly  developed  pro¬ 
cedure;  and  if  the  proper  input  values  art:  available  (e.g.,  line  param¬ 
eters,  laser  position,  molecular  concentrations),  it  generally  leads  to 
very  accurate  results.  The  procedure  is  one  of  simply  summing  the  con¬ 
tributions  to  the  LoLjI  absorption  coefficient  due  to  all  the  absorption 
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lines  i  for  each  molecular  species  j: 

m  n 

k  -  2  2  ku  <i2> 

J-l  i-l 

where  k_  is  the  absorption  coefficient  as  given  by  the  Lorentz,  Voigt, 
or  Doppler  expressions  (depending  on  the  total  pressure  as  discussed 
earlier  in  this  section;  for  the  ith  absorption  line  of  the  1th  molecular 
species.  Thus,  the  model  simply  evaluates  a  large  number  oi  individual 
absorption  coefficients  based  on  equations  given  earlier  and  a  line  pa¬ 
rameter  compilation  and  totals  them  all.  The  determination  of  how  many 
absorption  lines  n  to  consider  is  generally  based  on  the  location  of  the 
absorption  line  relative  to  that  of  the  laser.  That  is,  absorption  lines 
are  considered  one  by  one,  starting  at  the  nearest  to  the  laser,  until 
all  lines  within  seme  preselected  cutoff  distance  have  been  considered. 
This  cutoff  is  generally  taken  to  be  25  cm  ^  unless  there  is  information 
to  the  contrary.  Note  that  this  sumnatiori  is  carried  out  in  both  spec¬ 
tral  directions  from  the  laser  line. 

Several  computer  codes  have  been  developed  to  carry  out  this  cal¬ 
culation,  most  of  them  using  the  AFCL  Line  Parameter  Compilation  [3]. 
Probably  the  code  that  is  best  documented  and  most  widely  accepted  is 
the  AFCL  code  LASFR  {42],  This  code  Is  also  most  likely  to  be  maintained 
and  updated,  since  such  maintenance  apparently  is  one  of  the  charters  of 
the  Air  Force  Geophysics  Laboratory.  This  code  is  also  convenient  be¬ 
cause  it  incorporates  the  most  generally  accepted  models  for  molecular 
scattering  and  continuum  absorption,  as  well  as  typical  models  of  aero¬ 
sol  extinction. 

There  are  other  models  available  that  use  an  approximate  relation¬ 
ship  in  place  of  the  lorentz  format  large  distances  from  the  absorption 
line  center.  These  codes  (e.g.,  SYNSPF.C  of  Science  Applications,  Inc.) 
are  considerably  faster,  with  little-  degradation  in  accuracy;  so  con¬ 
sideration  of  absorption  line’s  far  from  the  laser  line,  as  necessary 
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in  some  spectral  regions  (e.g.,  CO,  HF  as  discussed  previously), 
can  be  accomplished  economically.  However,  SYNSPEC  is  neither  docu¬ 
mented  nor  widely  distributed,  so  it  would  probably  not  be  the  best 
choice  for  the  present  purposes. 

Another  class  of  model  has  been  developed  by  Tuer  [43]  which  is 
intended  for  easy  evaluation  of  molecular  absorption  for  a  limited  num¬ 
ber  of  laser  lines.  This  model  is  based  directly  on  the  results  of 
other,  more  sophisticated  computer  codes,  such  as  LASER  or  SYNSPEC, 
by  using  simple  analytic  functions  of  temperature  (T)  and  humidity 
fitted  to  their  calculated  results: 

k  1  ko + (v + v2)t + fa + vy  <i3> 

where  p  is  the  water  vapor  partial  pressure,  and  k^,  k^,  .  .  . ,  k^  are 
coefficients  determined  from  the  least-squares  fit.  An  example  of  the 
results  of  this  model  is  shown  in  Figure  8.  It  is  anticipated  that  such 
a  procedure  could  easily  and  effectively  be  applied  to  the  problem  of 
evaluating  the  effect  of  molecules  in  the  atmosphere  on  Jaser  transmis¬ 
sion,  for  specifying  safety  considerations.  Separate  tables  of  coef¬ 
ficients  would  be  required  for  various  total  pressures  so  that  the  ex¬ 
tinction  for  pa t>is  at  altitude  could  be  estimated. 

Molecular  Continuum  Absorption — In  certain  spectral  regions,  the 
primary  contributor  to  the  atmospheric  molecular  absorption  of  laser 
radiation  is  the  continuum  absorption  by  atmospheric  gases.  Therefore, 
the  accuracy  with  which  the  continuum  absorption  can  be  predicted  in 
these  spectral  regions  has  a  significant  impact  on  determining  the  safe 
laser  radiant  emittance  levels  for  a  given  situation.  The  accuracy  with 
which  the  infrared  transmission  in  the  atmosphere  can  be  calculated  has 
improved  slowly  since  the  development  in  the  early  70's  of  the  line-by¬ 
line  modeling  codes  by  McClatchey  et  al.  [2 1.  This  is  due  largely  to 
the  uncertainties  In  the  level  of  the  continuum  in  several  spectral  re¬ 
gions,  and  how  it  varies  with  temperature,  pressure,  and  molecular  con¬ 
centration. 
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Figure  8.  Sample  results  from  model  of  laser  radiation  as  a  function  of  temperature 
and  humidity,  for  Pi (9)  DF  line  at  2691.607  cm“* . 


In  this  section,  the  current  state-of -knowledge  on  theoretical  and 
experimental  aspects  of  infrared  molecular  extinction  is  presented.  Ex¬ 
periments  and  data  on  the  water  vapor  continuum  absorption  are  discussed 
also.  Included  is  a  brief  summary  of  as  yet  unpublished  findings  by 
Burch  in  the  3-5-pm  region  at  ambient  temperature.  Recent  nitrogen  and 
carbon  dioxide  continuum  absorption  experiments  are  discussed  below. 
Current  experiments  and  models  in  the  areas  of  continuum  absorption  by 
1^0,  ,  and  CC>2  are  also  compared.  Also  discussed  are  the  following 

continuum  absorption  mechanisms:  (1)  The  combined  effects  of  the  far 
wings  of  a  large  number  of  strong  lines;  (2)  transitions  within  dimers 
and  larger  polymers  of  water  vapor  molecules,  possibly  present  in  small 
concentrations  in  atmospheric  paths;  and  (3)  the  effects  of  far  wings, 
such  as  sub-Lorentzian  character  of  CC^  lines,  and  the  self -broadening 
effects  of  water  vapor  on  far  line  wings.  Finally,  in  Appendix  A,  re¬ 
cent  measurements  of  molecular  continuum  absorption  are  summarized. 

Temperature  dependence  of  the  continuum  absorption  of  water  vapor 
and  nitrogen  and  accurate  determination  of  the  spectral  shape  of  the  CC^ 
continuum  absorption  are  areas  of  major  interest  in  infrared  laser  ab¬ 
sorption  experiments.  Data  are  available  for  spectral  regions  3-5  urn 
and  8-12.5  pm  for  water  vapor.  3. 7-4. 8  pm  for  and  5-12.9  pm  and 
1.41-1.47  pm  for  CC^.  Of  primary  concern  in  current  experiments  is  the 
temperature  dependence  of  the  (^0  continuum  absorption  in  the  3-5-pm  and 
8-12.5— Pm  regions.  Reliable  quantitative  data  in  these  spectral  regions 
at  realistic  atmospheric  conditions  will  lead  to  more  reliable  attenua¬ 
tion  predictions,  and  possibly  to  an  understanding  of  the  mechanism  re¬ 
sponsible  for  the  absorption  in  these  regions. 

The  continuum  absorption  in  the  3-5-pm  spectral  region  is  thought 
to  be  simply  far-wing  absorption  of  the  strong  water  vapoj;.  lines  that 
are  a  few  hundred  wavenumbers  away  [37],  The  investigation  of  the  pres¬ 
sure  and  temperature  variation  of  the  weak  vapor  absorption  in  the  3-5- 
pm  window  has  tended  to  follow  the  experimental  approach  of  Burch  et  al. 
[44],  Cosden  et  al.  [37]  present  data  in  the  3-5-pm  region  on  the 
temperature  and  pressure  variation  of  the  self-absorption  coefficient  of 


H^O,  and  of  tin.”  ratio  of  the  ioreign  to  self-broadening  coefficients  at 

three  temperatures  —  338,  384,  and  428  K.  Briefly,  the  conclusions  of 

their  work  were:  (1)  Absorption  due  to  water  vapor  in  the  far  wings  is 

2 

quadratic  in  the  water  vapor  pressure  (i.e.,  k  a  p  );  (2)  In  the  presence 
of  a  foreign  gas,  additional  "broadening"  is  observed  which  varies  lin¬ 
early  with  the  foreign  gas  pressure,  giving  the  continuum  absorption  co¬ 
efficient  in  this  region  the  form: 


where 

pl}  Q  =  Pressure  in  atmosphere  of  water  vapor 

=  Pressure  in  atmosphere  of  water  vapor  and  foreign  gas 

C  =  Wavelength-dependent  coefficient  obtained  by  Burch 
2  (see  Figure  9) 

B  =  Foreign  broadening  coefficient; 

and  (3)  Continuum  absorption  decreases  as  the  temperature  (T)  is  in¬ 
creased,  following  the  form  exp(-const./T) . 

Since  this  early  work,  major  efforts  have  been  made  by  Damon  et 
al.  145],  Mills  120],  White  et  al,  [46,47],  and  Watkins  et  al.  [10] 
to  verify  and  extend  Burch's  data.  These  groups  utilized  deuterium 
fluoride  (OF)  lasers  having  26  lines  in  the  3-5-fim  window.  Damon  et  al. 
[45]  present  spectrophone  measurements  of  the  absorption  coefficient  as 
a  function  of  water  vapor  partial  pressure  up  to  approximately  15  torr, 
with  each  sample  buffered  to  a  total  pressure  of  760  torr  with  artificial 
air  (80  percent  N^/20  percent  0^).  To  extract  H^O  continuum  absorp¬ 
tion  coefficient  at  14.3  torr,  the  measured  11D0  absorption  and  the  cal¬ 
culated  (AFGL  tape)  H.;0  absorption  from  nearby  lines  were  subtracted 
from  the  total  extinction.  The  results  for  the  absorption  coefficient 
at  six  wavelengths  are  reproduced  in  Table  10  and  plotted  in  Figure  10. 
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Figure  10.  Measured  H2O  continuum  absorption  coefficients 
as  a  function  of  frequency  compared  to  Burch 
extrapolation. 


These  authors  quote  a  water  vapor  continuum  absorption  coefficient  at 
these  wavelengths  of  k»0.04i0.02  cm  however,  the  data  indicate  a 
somewhat  larger  spread  in  uncertainty.  Note  that  this  is  for  isotropic¬ 
ally  pure  1^0  at  a  partial  pressure  of  14.3  torr  and  a  temperature  of 
297  K.  In  Figure  10,  we  have  compared  their  results  with  an  absorption 
coefficient  obtained  by  extrapolating  the  Burch  data  given  in  Figure  9. 
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Using  the  values  of  the  self -broadening  coefficient  from  Figure  9  in 
Eq.  14,  we  obtain  extinction  coefficients  that  range  between  0.02  and 
0.03  km  ^  in  the  wavenumber  range  of  2600  to  2800  cm  .  The  agreement 
in  this  spectral  range  is  good  and  lends  support  to  the  Burch  extrapo¬ 
lation;  however,  the  support  is  weak  at  best  because  the  quoted  error  in 
the  spectrophone  data  is  so  large  (±50  percent). 

Also  shown  in  Table  10  and  Figure  10  are  data  from  two  White  cell 
(1.34-km  pathlength)  measurements,  by  Damon  [45]  and  Mills  [20],  at 
14.3  torr  water  vapor  pressure,  760  torr  total  pressure  (buffered  with 
,  and  at  a  temperature  of  297  K.  The  uncertainty  presented  in  Table 
10  for  the  spectrophone  data  is  that  quoted  by  the  experimenters  for 
their  average  absorption  coefficient.  Mill's  White  cell  and  Damon's 
spectrophone  values  differ  from  each  other  by  2.5  standard  deviations 
at  the  lower  frequency.  At  the  higher  frequency,  the  two  measurements 
are  in  better  agreement,  although  they  are  not  within  one  standard  devi¬ 
ation. 


The  most  recent  data  available  on  water  vapor  absorption  in  the 

3-5  pm  window  are  from  White  et  al.  [46,47]  and  Watkins  et  al.  [10]. 

White  et  al.  measured  the  water  vapor  absorption  coefficient  at  25  DF 

lines  in  a  White  cell  at  338  K,  with  a  water  vapor  partial  pressure  of 

72  torr.  White  cell  and  spectrophone  data  are  also  presented  for  T»296  K 

and  14.3  torr  water  vapor.  All  samples  were  buffered  to  a  total  pressure 

of  760  torr  with  a  4:1  mixture  of  ^  to  O2.  The  nitrogen  continuum  was 

experimentally  subtracted.  The  338  K,  ptI  »  72  torr  data  were  taken  to 

rl2^ 

allow  a  direct  and  unambiguous  comparison  with  Burch's  data  taken  under 
similar  conditions,  and  to  provide  data  beyond  the  frequency  range  cov¬ 
ered  by  Burch  at  this  temperature  and  pressure.  The  continuum  absorption 
coefficients  were  derived  by  subtracting  the  HDO  and  1^0  line  contribu¬ 
tions  from  the  total  water  vapor  absorption.  Their  high-temperature  re¬ 
sults  are  compared  in  Figure  11.  The  calculated  HDO  and  H2O  line  con¬ 
tributions  were  based  on  the  January  1976  updated  AFGL  data  tape.  Con¬ 
sidering  the  30  percent  uncertainty  contained  in  the  Burch  results,  the 
two  sets  of  data  in  Figure  11  are  consistent  in  this  high-temperature, 
high-pressure  region. 


42 


TABLE  10.  H2O  EXTINCTION  COEFFICIENT  MEASURED  BY  DAMON  [45]  AND  MILLS  [20] 
WITH  CALCULATED  H2O  LINE  ABSORPTION  AND  MEASURED  HDO  ABSORPTION 
SUBTRACTED  (pu  -  -  14.3  torr;  p_  -  760  torr;  T  -  297  K) 


0.2S 


FREQUENCY  fan"1) 

Figure  11.  Water  vapor  continuum  measured  by  White  et  al. 

[46,47]  and  Burch  [44]. 

Additional  continuum  data  by  White  et  al.  [47]  are  presented  and 
compared  in  Figure  12  with  the  Burch  extrapolations  for  nominal  atmo¬ 
spheric  conditions  (i.e.,  'f-296  K  and  p..  -14.3  torr) .  One  set,  from 
spectrophone  and  White  cell  measurements,  represents  absorption  by 
natural  water  vapor,  from  which  the  AFGL  calculated  !UX)  and  i^O  line 
contribution  have  been  subtracted.  The  second  set  is  White  cell  data 
from  HDO-depleted  (2  percent  of  natural  abundance)  water  vapor  samples. 
For  ease  in  comparing  data,  the  error  bars  on.  the  HM)-depleted  data  and 
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•CALCULATED  HDO  LINE  ABSORPTION  HAS  BEEN  SUBTRACTED. 

Figure  12.  Water  vapor  continuum  measured  by  White  et  el.  [46,47] 
compared  with  the  extrapolated  Burch  continuum; 

T  -  296  K,  p„  .  ■  14.3  torr. 
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the  spectrophone  data  were  omitted.  Their  error  bars  were  comparable  to 
those  shown.  These  data  were  averaged  by  White  et  ai.  and  then  com¬ 
pared  to  the  Burch  result.  From  this  comparison,  they  concluded  that 
the  Burch  extrapolation  underestimates  the  continuum  in  the  3-5-pm  region. 

We  conclude  the  discussion  of  the  state-of-knowledge  in  3-5-pm 
experiments  with  a  summary  of  the  sequel  to  this  last  work  (by  White) , 
in  which  Watkins  et  al.  [  1 0]  present  new  White  cell  data  at  298  K, 

14.3  torr  water  vapor,  and  several  different  foreign-broadening  gac 
pressures.  Their  data  for  764.3  torr  total  pressure  is  in  general 
agreement  with  the  earlier  data  by  White  [46]  (i.e.,  falls  within  the 
error  bars  of  White's  earlier  data  presented  in  Figure  12),  but  lies 
systematically  closer  to  the  Burch  extrapolation  at  this  temperature. 

These  authors  also  present  (see  Figure  13)  new  determinations  for  the 
ratio  of  foreign-  to  self-broadening  coefficients  (Cf/Cg)  at  26  DF  laser 
frequencies.  Their  weighted  average  value  for  this  ratio  is  C^/Cg=0.011, 
at  1*298  K  and  pH^0=14.3  torr.  The  weighted  average  is  formed  by  weight¬ 
ing  the  individually  measured  values  of  C^/Cg  with  the  measured  extinc¬ 
tion  coefficients  at  zero -broadening  pressure  and  normalizing  to  the 
sum  of  extinction  coefficients.  This  weighted  average  is  expressed  in 
the  following  form: 


where  the  k(p^*0)  is  the  extinction  coefficient  extrapolated  to  zero 
foreign  gas  pressure.  This  ratio  of  foreign-  to  self-broadening  is  more 
than  ..a  order  of  magnitude  smaller  than  the  value  obtained  by  Burch  [44] 
at  338  K  (i.e.,  Cf/Cs»0.12) . 

Based  on  their  result  for  the  ratio  C  /C  ,  Watkins  et  al.  propose 
a  somewhat  ad  he<"  model  for  the  water  vapor  continuum  absorption  coef¬ 
ficient  at  these  frequencies.  By  assuming  a  third  contribution  to  the 
continuum  absorption  coefficient,  which  is  independent  of  foreign- 
broade.'ng  pressure,  and  reanalyzing  their  data,  they  obtained  the  value 
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Figure  13.  Ratio  of  foreign-  to  self-broadening  coefficients 
at  26  DF  laser  lines  before  (•)  and  after  (a) 
curve  fitting  the  measured  absorption  coef t ic.ants, 
Watkins  et  al  .  [10], 


of  this  third  contribution.  This  additional  absorption  term  is  hypoth¬ 
esized  to  result  from  water  vapor  dimers.  With  this  assumption  and 
the  additional  assumption  that  the  Burch  extrapolation  is  valid  and 
represents  only  far-wing  absorption,  the  Burch  continuum  is  subtracted 
from  their  data  at  each  wavelength.  So  by  hypothesis,  wnat  remains  is 
the  water  vapor  dimer  absorption  coefficien*.  The  temperature  depen¬ 
dence  of  this  tern  is  under  investigation  by  Watkins  et  al.  [10], 

Burch  has  completed  a  new  experiment  in  the  3-5-pm  region  on  the 
water  vapor  continuum  at  both  high  and  ambient  temperatures  [48].  His 
high-temperature  data  (338  K)  agree  with  his  earlier  data,  and  hence 
with  the  338-K  data  of  White  et  al.  [47 j.  His  low-temperature  (near 
296  K)  data,  however,  fall  below  those  of  White  et  al.  Burch  has 
found  that  heating  the  mirrors  in  White  cells  during  ambient  temperature 
experiments  can  produce  increases  of  up  to  an  order  of  magnitude  in  the 
transmittance,  lie  attributes  this  to  absorbed  water  vapor  on  the  mirror 
surfaces.  McClatchey  indicated  [48]  that  the  new  version  of  LOWTRAN, 
LOWTRAN  V,  will  use  a  water  vapor  continuum  in  the  3-5-pm  region  based 
on  reconmendations  by  Burch.  The  new  model  will  probably  be  lower  than 
White's  data  [49]. 

In  concluding  the  state  of  the  3-5-pm  water  vapor  continuum  mea¬ 
surements,  it  is  worth  noting  that  in  1971,  Burch  [44]  reported  that 
the  transmitted  signal  slowly  decreased  over  a  period  of  several  hours 
following  the  filling  of  his  White  cell.  He  attril>uted  at  least  part 
of  this  anomalous  absorption  to  a  water  vapor  film  that  formed  on  the 
mirror  surfaces.  He  also  found  that  the  signal  could  be  retrieved  by 
complete  evacuation  of  the  cell.  These  early  results  of  Burch  and  his 
recent  attention  to  the  temperature  dependence  of  this  absorption  lead 
to  the  reasonable  conclusion  that  until  this  large  systematic  error  is 
quantified,  good  low-temperature  White  cell  data  on  water  vapor  contin¬ 
uum  absorption  are  accurate  to  approximately  *30  percent  in  this  wave¬ 
length  region.  This  is  the  error  quoted  by  Burch  [44]  in  1971  for  his 
low-temperature  data. 


48 


In  the  8-12-pm  region,  early  field  data  by  Kflndrat'yev  et  al.  [50], 
taken  at  three  elevations  (sea  level,  310  m,  and  3100  m)  under  carefully 
monitored  atmospnei ic  conditions,  show  the  general  dependence  of  the  con¬ 
tinuum  absorption  coefficient  on  wavelength  (Figure  14).  Considering  the 
lack  of  control  over  such  factors  as  temperature  and  aerosol  content  in¬ 
herent  in  their  experiment,  the  data  can  only  be  regarded  as  qualitative. 

One  of  the  earliest  experiments  claiming  to  provide  evidence  for 
the  H?0  dimer  contribution  to  continuum  absorption  in  the  8-12-um  window 
is  that  of  Varanisi  [51].  The  absorption  coefficient  for  PH2q=:'  atK  at 
three  rather  high  temperatures  (400,  450,  and  500  If)  is  plotted  in  Figure 
15  as  a  function  of  frequency.  To  obtain  a  semi  empirical  estimate  of  the 
dimer-binding  energy,  Varanisi  first  defines  what  he  terms  "an  average 
absorption  coefficient  per  particle"  by  the  relation: 


WAVELENGTH  (gml 


Figure  14.  i)ual  Itaf  lve  sliapc  of  wavelength  dependence  of  H2® 
continuum  abso,ptlon  coefficient  from  data  of 
Kondrat’yev  et  al .  [5()J. 
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Figure  15.  Plot  of  the  absorption  coefficient  from  Varanisi 

[51]  for  p„  “2  atm  at  three  elevated  temperatures. 

where  R  is  the  gas  constant  per  gram.  Using  tabula  d  values  (Benedict 
and  Kaplan  [52])  for  the  line  intensities  with  a  temperature  dependence 
assumed  to  be  described  by: 


o  a  N  -3/2  -(E/RT) 

S  “  -  T  e 


-3/2 

where  T  is  the  temperature  dependence  of  the  rotational  partition 

3/2  _  _i 

function,  Varanisi  plots  £n(T  ['^Theory}  versus  T  •  The  tabulated 
line  intensities  used  are  for  temperatures  of  220,  260,  and  300  K.  Thi 


(ft* 


plot  yields  a  straight  line  which  can  be  extrapolated  to  obtain  a  theo¬ 
retical  estimate  <>f  .n  average  absorpt ion  coefficient,  {kTj,_,  ,  at  the 

experimental  temperatures  of  Interest.  We  assume  that  the  difference 

fi(tT>  r  (lT)„easu,ed  -  '^Theory  (18> 

between  the  measured  average  absorption  coefficient  (reduction  of  this 
quantity  from  Varanisi’s  data  is  unclear)  and  the  theoretical  average 
absorption  coef  t  re i*_r.t  is  proportional  to  the  number  of  hydrogen  bonds 
with  a  Boltzmann  temperature  dependence: 


U  kT)u 


-E„-H/RT 


(19) 


Then  a  value  of  5  kcal/mole  is  derived  for  the  h/drogen-bonding  energy, 
E„  from  the  data  plotted  in  Figure  15  (taken  at  p„  =  2  atm).  Bv  in- 

rt— n  n2^ 

creasing  the  pressure  to  pu  ,=10  atm,  a  value  of  3  kcal/mole  was  deter- 
mined.  Varanist  concludes  by  asserting  tliat  his  results  for  do  not 

change  significantly  when  allowance  for  far- wing  absorption  of  strong 
distant  lines  is  made. 


McCoy  et  al.  13.  present  White  cell  data  on  foreign-  and  self¬ 
broadening  H^O  absorption  at  10.59  and  9.55  pm  at  a  temperature  of  298  K. 

A  least-squares  fit  at  a-10.59  um  to  a  plot  of  -P,n(transmittance)  versus 
2 

p  .  was  used  to  obtain  a  self-broadened  extinction  coefficient,  given 

'2°  2 
as  a  function  of  p^  Q  by 


k( 10.59  urn)  -  8.39xlO~4  pu  ^(knf1)  . 

H2° 


(20) 


Their  estimated  error  was  *0.587x10  In  the  coefficient;  p  was  in 

1*2  U 

torr.  At  9.55  um,  these  authors  obtained 


k(9. 55  pm)  -  b.n7xlO~4  pjj  0(km-1) 


(21) 
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with  an  error  of  10.8x10  .  Further  data  are  presented  for  the  10.59-pm 

transmittance  in  the  presence  of  foreign  gas.  In  these  experiments,  the 
H.,0  vapc-r  pressure  was  varied  while  a  constant  total  pressure  of  700  torr 
was  maintained.  A  fit  of  the  equation 

k(10. 59  um)  =  opH  0Jpt+(B-1)ph  Qj  (22) 


to  their  transmittance  data  yielded 
a  =  4 . 32xl0~6 

B  =  Self-  to  foreign-broadening  coefficient  =  194 
=  Total  pressure 

with  an  estimated  error  of  HO  percent  in  a.  L'sing  these  coefficients 
and  expressing  the  attenuation  in  units  of  dB-km  at  a  total  pressure 
of  760  torr,  the  absorption  loss  at  sea  level  was  given  by  McCoy  et  al. 
as 


Loss  =  1.43xl0~2  pu  _  +  3.62x10  3  p2  dE-km  1  .  (23) 

t^O  H20 

To  test  this  result,  field  data  were  taken  at  10.59  pm  over  a  1.95-km 
outdoor  path.  The  measurements  were  made  at  temperatures  ranging  from 
15°  to  30°C,  and  relative  humidities  of  50  to  90  percent.  The  field 
data  are  plotted  in  Figure  16,  along  with  the  calculated  loss  using 
Eq.  23  plus  a  constant  0.47  dB-km  1  (445  ppm  CO 2)  loss  due  to  C02  ab¬ 
sorption  (solid  curve) .  The  authors  point  out  that  the  445  ppm  C02  is 
higher  than  the  330  ppm  C02  (see  dashed  line  in  Figure  16)  usually  as¬ 
sumed  for  a  standard  25°C  atmosphere;  however,  this  higher-than-average 
C02  content  was  substantiated  by  measurement  of  the  C02  concentration. 

By  vircue  of  the  agreement  between  the  field  data  and  the  calculated 
attenuation  based  on  Eq.  23,  McCoy  et  al.  conclude  that  the  absorption 
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Figure  16.  Field  data  of  McCoy,  et  al.  [13],  at  10.59  pm  com¬ 
pared  with  calculated  attenuation  based  on  Eq.  23 
text. 

loss  is  not  highly  temperature  sensitive  over  the  temperature  range  of 
their  data  (15°  to  30°C). 

Bignell  [53]  presents  water  vapor  absorption  at  several  wave- 
lengths  in  the  8-12-pm  range.  To  extract  the  quoted  values  of  the  self¬ 
broadening  coefficient,  the  author  removed  the  foreign-broadening  con¬ 
tribution  from  his  raw  data  without  quoting  the  values  or  the  source  of 
the  foreign-broadening  coefficient  used. 
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Measured  absorption  coefficients  by  Moskalenko  et  al.  [54],  at 
10.6  Mm  as  a  function  of  temperature  are  displayed  in  Figure  17.  These 
authors  give  no  interpretation  or  possible  explanation  for  their  obser¬ 
vations.  The  behavior  of  the  water  vapor  absorption  coefficient  as  a 
function  of  temperature  reported  by  these  authors  has  not  to  date  been 
corroborated  by  other  researchers  in  this  field. 

In  1976,  Gryvnak  et  al.  [55]  reported  data  at  three  temperatures 

in  the  8-12-Mm  region,  which  modified  slightly  their  earlier  results 

[39]  for  C°,  the  water  vapor  self -broadening  coefficient.  In  some  cases, 

the  changes  in  C°  were  as  great  as  20  percent.  These  later  results  are 
s 

plotted  in  Figure  18  for  all  three  temperatures  —  296,  392,  and  430  K. 

The  authors  quote  an  upper  bound  of  ±10  percent  for  the  error  in  the  392 

and  430  K  data.  The  sharp  increase  in  C°  with  increasing  frequency  above 
-X  s 

1150  cm  (8.70  Mm)  is  attributed  by  these  authors  to  an  increasing  con¬ 
tribution  from  the  lines  centered  just  above  this  frequency.  However, 


Figure  17.  Measured  water  vapor  absorption  coefficients  as 
a  function  of  temperature  [54], 
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FREQUENCY  (on-1) 

Figure  18.  Water  vapor  self-broadening  coefficient,  C®,  as 
a  function  of  wavelength  at  three  temperatures 
[55]. 

they  assert  that  a  large  portion  of  the  continuum  below  1100  cm  ^  (above 

9.1  Um)  is  probably  due  to  the  extreme  wings  of  the  very  strong  lines 

centered  below  600  cm  ^  (above  16.7  pm).  This  could  account  for  the 

increase  in  the  self-broadening  contribution  with  decreasing  wavenumber 

below  1100  cm  The  conclusion  was  unchanged  from  the  earlier  (1970) 

report  by  these  authors  [39]  that  the  ratio  of  foreign-  to  self-broadening 

coefficient,  C°/C°,  is  less  than  0.005  near  296  K, 

’  n  s’ 

Shumate  et  al.  [56]  measured  extinction  coefficients  at  49 
wavelengths  within  two  laser  bands,  one  centered  at  9.4  pm  and  one 
at  10.4  pm.  These  measurements  were  taken  with  a  spectrophone  at  three 
water  vapor  partial  pressures  buffered  with  air  to  a  total  pressure  of 
76C  torr,  at  a  temperature  of  300  K  (see  Tables  11  and  12).  Direct  com¬ 
parison  with  the  measurements  of  McCoy  et  al.  [13]  at  10.59  pr.i,  dis¬ 
cussed  earlier,  is  given  m  Table  13.  The  agreement  is  fair  according 
to  Shumate  et  al.  We  can  compare  the  data  of  Gryvnak  et  al .  [55]  in 
this  wavelength  region  with  the  data  from  Shumate  et  al.  To  make  a 


TABLE  11.  ABSORPTION  COEFFICIENTS  OF  WATER  VAPOR  IN  AIR 
AT  THREE  WATER  VAPOR  PARTIAL  PRESSURES: 
10.4-ym  BAND  OF  THE  CX20j'*  LASER  [56] 


Laser 

line 

(001-1) 

Frequency 

(cm-) 

A.Iuorption  coefficients*  (10** 

cm—) 

5.0  Tore 

10.0  Tore 

1 5.0  Tonr 

P(  32) 

932.960 

'1.39,*  0.30* 

2.32,*  0.82* 

3.63*  1.93* 

A  ( 30) 

934.898 

0.36,*  0.30* 

0.95.*  0.87* 

2.01,*  1.92* 

P(28) 

936.804 

0.33,  0.30* 

0.87 

1.85 

P(  26) 

938.688 

0.37,  0.34* 

0.89 

1.78 

*(24) 

940.648 

0.36 

1.00 

1.92 

A0!  22) 

942.383 

0.37 

0.94 

1.97 

1»(20) 

944.194 

0.41 

1.15 

2.18 

PU8) 

948.980 

0.38 

1.00 

1.97 

*(16) 

947.742 

0.49, *■<*  0.45*.* 

.1.23,*  1.18* 

2.66.*  2.60* 

P(14) 

949.479 

0.45,*  0.40* 

1.27*  1.19* 

2.50*  2.41* 

*(12) 

981.192 

0.35*  0.31* 

1.01.*  0.95* 

2.10.*  2.03* 

*(10) 

952.881 

0.32 

0.91 

2.19 

*(8) 

967.707 

1.96*  0.41* 

3.08*  0.88* 

4.53.*  2.03* 

*(10) 

969.139 

0.40,*  0.34* 

0.81.*  0.74* 

1.55,*  1.46* 

*(12) 

970.847 

1.00,*.*  0.97*.* 

2.12,*  2.07* 

3.32.*  3.27* 

*(14) 

971.930 

0.97  *.*  0.92'  -* 

1.81,*  1.75* 

2.63,*  2.56* 

*(  16) 

973.289 

0.78* 

1.70 

2.59 

*(18) 

974  622 

0.61* 

1.22 

1.88 

71(20) 

975.930 

5.49* 

11.0 

16.8 

*(22) 

977.214 

0.84* 

1.66 

2.67 

*(24) 

978.472 

0.50 

1.12 

1.90 

*(  26) 

979.705 

0.49 

1.03 

1.78 

*(28) 

980.913 

0.40 

,  0.92 

1.66 

*(30) 

982.096 

0.46* 

1.07 

1.97 

71(32) 

983.252 

0.40 

1.09 

1.96 

"Measured 

at  760-torr  total  pressure 

,  300 -K  temperature. 

^Significant  absorption  from  ammonia  observed  at  this  frequency. 

cCorrected  by  subtracting  the  ammonia  contribution. 

^Significant  absorption  from  water  line?  predicted  at  this 
frequency. 


TABLE  12.  ABSORPTION  COEFFICIENTS  OF  WATER  VAPOR  IN  AIR 
AT  THREE  WATER  VAPOR  PARTIAL  PRESSURES: 

9.4-yto  BAND  OF  THE  C12o£‘  LASER  [561 


Liter 

ft.  ,  A  n  .  1  A  tt  e  tl 

Abtorption  coefficient**  (10“* 

cm'*) 

lint 

(001  -11 ) 

r lequrncy 

(cm-) 

5.0  Torr 

10.0  Tore 

15.0  Tore 

P<32) 

1035  474 

0.42.*  0  39* 

0  85.*  0.82* 

1-47.*  1.43* 

P(30) 

1037.434 

0.39 

0.77 

1.39 

P<28) 

1039.369 

1.30" 

2.72" 

4.46" 

P{26) 

1041.279 

0  41 

0.77 

1.37 

P(24) 

1043 163 

0.41,*  0.37* 

0.93,*  0  89* 

1.57,*  1.52* 

P(22) 

1045.022 

0.47,*  0.44* 

0  96,*  0.92* 

1.62,*  1.58* 

P(20) 

1046 854 

0.51,*  0.36* 

1.09,*  0.88* 

1.69,*  1.44* 

P(18) 

1048  661 

0.37,*  0  35* 

0  89,*  0.86* 

1.43,*  1.40* 

/*(16) 

1050.441 

0.47, »•"  0  45*  " 

1.02, *-"1.00*." 

1.58,*"  1.56*" 

P(U) 

1052.196 

0.41,*  0.38* 

1.03,*  1.00* 

1.58.*  1.55* 

P(12) 

1053  924 

0  43,*  0.36* 

0  86.*  0.77* 

1.50,*  1.40* 

P(10) 

1055.625 

1.15" 

2.82" 

3.86" 

R(10) 

1071  884 

0.35 

0.83 

1.25 

F(12) 

1073  278 

0  44" 

1.17" 

1.88" 

fHU) 

1074  646 

1.48" 

3.51" 

5.24" 

R(  16) 

1075  988 

1.53.*-"  0.57*." 

2. 88,*"  1.60*" 

3.54.*"  2.10*" 

fi(18) 

1077.303 

0  43 

0.99 

1.45 

R(20) 

1078  591 

0  40 

0  93 

1.41 

R(  22) 

1079  852 

0  40 

0.91 

1.38 

R(24) 

1081  087 

0  39 

0.95 

1.47 

R(26) 

1082  296 

0  36 

0.84 

1.35 

R(28) 

1083  479 

0.41,*  0.38* 

0.97,*  0.93* 

1.42.*  1  38* 

Rl  30) 

1084  635 

6.21.*  1.40* 

8  25,*  1.65* 

9.36.*  1.93* 

R(32) 

1085  765 

0  58" 

1.27" 

1.85" 

Measured  at  760-torr  total  pressure,  300-K  temperature. 

^Significant  absorption  from  ammonia  observed  at  this  frequency. 

£ 

Corrected  by  subtracting  the  ammonia  contribution. 

^Significant  absorption  from  water  lines  predicted  at  this 
frequency. 


TABLE  13.  COMPARISON  OF  MEASURED  WATER  VaPOR  ABSORPTION 
COEFFICIENTS  AT  10.59  pm 

p  Absorption  coefficient  (km 

h2o  - 


(torr) 

Shumate  et  al.  [56] 

McCoy  et  al.  [13] 

5 

0.041 

0.038 

10 

0.115 

0.117 

15 

0.218 

0.235 

20 

0.353 

0.396 

comparison,  we  assume  C°/C“=0.005  [55],  and  p^  ^=15  torr  with  a  total 
pressure  of  760  torr.  This  comparison  is  given  in  Table  14.  While  the 
agreement  is  poor,  it  must  be  remembered  that  the  comparison  is  dependent 
upon  the  assumption  of  a  constant  foreign-  to  self-broadening  ratio  of..  ^ 
0.005.  A  choice  of  smaller  values  would  naturally  yield  a  more  favorable 
comparison. 


TABLE  14.  COMPARISON  OF  MEASURED  WATER  VAPOR  ABSORPTION 
COEFFICIENTS  AROUND  10  pm  (Ph,o’15  torr* 

P  l«760  torr,  T-300  K) 


Frequency 

Absorption  coefficient 

(km  ) 

(cm-1) 

Gryvnak  et  al.  [55]  Shumate 

et  al.  [56] 

935 

0.274 

0.192 

981 

0.262 

0.166 

1072 

0.190 

0.125 

Coffey  [57]  presents  fie Id -measured  water  vapor  attenuation  at 
10.7  and  11.6  pm,  recorded  during  flights  over  sea  environments.  He 
analyzed  his  results  by  writing  the  water  vapor  absorption  coefficient 
in  terms  of  the  observed  transmittance  t  as  follows: 

k(g  ^cm^)  •  — 1/u  2n  T  (24) 
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where  t  is  the  measured  transmission  and  u  is  the  absorber  amount. 
Coffey  then  writes  k  as  the  sum  of  two  terms: 


k  =  k 


l(%  +  8Ph2o)  +  »H2Ok2 


where  p  and  Ph2o  are  t*ie  nitrogen  and  water  vapor  partial  pressures  in 
atmospheres,  respectively,  and  3  is  the  self -broadening  coefficient,  as¬ 
sumed  by  Coffey  to  be  6,  based  on  Reference  58.  In  order  to  extract 
from  k  the  values  of  k2,  Coffey  performs  a  line-by-line  calculation  to 
obtain  k^  which  was  then  subtracted  from  the  data.  The  line-by-line  cal¬ 
culations  were  based  on  McClatchey's  1973  line  data  [3]  (portions  of 
these  HDO  line  parameters  underwent  major  revisions  in  1976).  The  gen¬ 
eral  trend  of  the  results  of  this  analysis  is  to  confirm  the  negative 
temperature  dependence  of  the  water  vapor  absorption  observed  by  Burch 
et  al.  [39].  Any  quantitative  comparison  would  suffer  from  uncertain¬ 
ties  in  (1)  the  computed  values  for  k^,  and  (2)  the  assumed  value  for 
the  ratio  of  foreign-  to  self-broadening  coefficients.  Comparing  the 
slope  from  a  £n(k2)  versus  T  ^  plot  with  the  temperature  coefficient 
obtained  from  the  second  virial  coefficient  of  steam,  Coffey  concludes, 
that  water  vapor  dimers  are  responsible  for  absorption  in  the  10-12~ym 
region. 

Montgomery  [59]  presents  measurements  of  the  temperature  depen¬ 
dence  of  the  self-broadening  coefficient,  Cg,  near  8.33  ym,  using  a 
lead-tin-telluride  diode  laser  and  a  40. 5  -  m  -pathlength  White  cell. 

The  temperature  used  ranged  from  333  to  473  K.  Montgomery  first  mea¬ 
sured  the  total  extinction  as  a  function  of  the  water  vapor  partial 
pressure  at  each  temperature.  A  typical  example  of  this  data  is  plotted 
in  Figure  19,  where  the  error  bars  indicate  fluctuations  in  laser  in¬ 
tensity  during  the  measurements.  The  straight  line  is  a  fit  of  the  data 

O 

t°  k/p^Q=Const.  The  result  of  this  fitting  at  each  temperature  is  then 
used  in  the  expression 

k/»H20  ■  Cs/kBolt.T  «« 
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WATER-VAPOR  PARTIAL 
PRESSURE,  PH2q  (atm) 


Figure  19.  Pressure  dependence  of  the  water  vapor  absorption 
coefficient  at  423  K,  from  Montgomery  [59], 


to  arrive  at  a  self-broadening  coefficient.  Montgomery’s  results  for 
Cg  as  a  function  of  temperature  are  plotted  in  Figure  20,  along  with 
values  from  Gryvnak  et  al.  [55]  for  comparison.  It  is  reasonable  to 
say  the  two  sets  of  data  are  in  general  agreement. 

Nordstrom  et  al.  [60]  present  water  vapor  absorption  measurements 
at  room  temperature  (295,5  ±  0.5  K)  for  five  CO^  laser  wavelengths  in 
the  10.4-ym  band.  Three  gases,  pure  nitrogen  and  two  nitrogen-oxygen 
mixtures  (80:20  and  60:40),  were  used  to  buffer  the  water  vapor  samples 
to  a  total  pressure  of  760  torr.  Transmission  measurements  were  made 
over  a  range  of  water  vapor  partial  pressures  from  0  to  15  torr  for  each 
buffer  gas.  These  authors  analyze  their  data  at  each  of  the  five  wave¬ 
lengths  in  terms  of  an  extinction  coefficient  of  the  form: 

k  "  W(CsPh20  +  Cnpn  +  Vo)  (27) 
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Figure  20.  Comparison  of  measured  water  vapor  self-broadening 
coefficients  near  8.33  ym  as  a  function  of  temper¬ 
ature. 


where  Cg  and  are  the  usual  self-  and  foreign-broadening  coefficients 
and  Cq  is  an  oxygen  foreign-broadening  coefficient.  Since  one  of  the 
laser  lines  [R(20),  A-io.2466  um}  is  so  very  near  a  water  line,  the 
authors  used  this  line  to  determine  Cq.  This  value  for  Cq  was  assumed 
to  apply  to  all  lines  investigated.  A  Lorentz  line  shape  with  a  width 
given  by: 

a  -  <*0[P+(C0/C-l)p|J]  (28) 
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where  Pq  is  the  average  oxygen  pressure  and  P  is  the  total  pressure, 
leads  to  an  absorption  coefficient  given  by: 


k(v  2  Vq) 


ph2o 


(29) 


A  least-squares  fit  to  k  versus  p^  yielded  a  value  for  the  ratio  Cq /C^ 
of  0.75;  i.e.. 


CjC  =  0.75 
0  n 


(30) 


which  Nordstrom  et  al.  interpret  as  strong  evidence  that  oxygen  is  not 
so  efficient  a  broadening  gas  as  nitrogen.  The  results  for  the  self¬ 
broadening  coefficient  and  the  ratios  of  the  foreign-  to  self-broadening 
coefficients  are  tabulated  in  Table  15  along  with  the  results  from 
Gryvnak  et  al.  [55 j  for  comparison.  These  authors  conclude  by  pointing 

out  that  the  small  values  of  C  /C  force  a  fit  of  Lorentzian  wings  to 

n  s 

the  observed  absorption;  i.e.,  the  Lorentz  line  shape  is  not  adequate 
to  model  absorption  in  the  wings  of  strong  water  lines. 


TABLE  15.  COMPARISON  OF  MEASURED  SELF-BROADENING 
COEiFICIENTS  AND  FOREIGN-  TO  SELF¬ 
BROADENING  RATIOS  FROM  NORDSTROM  et  al. 
[60]  WITH  THOSE  FROM  GRYVNAK  et  al.[55] 


V 

C* (cm2 -atm-1  - 
*  x  10” 

molecule-1) 

cVc* 

n  e 

******* 

Nordstrom 

Gryvnak 

Nordstrom 

Gryvnak 

936.804 

2.20 

1.25 

0.005 

0.005 

944.194 

2.11 

2.19 

0.003 

0.005 

977.214 

1.84 

2.18 

0.009 

0.005 

980.913 

1.75 

2.00 

0.003 

0.005 

The  most  recent  data  on  vater  vapor  nitrogen  attenuation  coef¬ 
ficients  are  from  Peterson  et  al .  [61  j,  who  present  self-broadening  co¬ 
efficients  and  ratios  C  /C  at  several  C0„  laser  frequencies.  Spectro- 

ns  i. 

phone  and  White  cell  data  were  concluded  to  be  in  good  agreement.  As 
the  spectrophone  data  were  not  at  a  single  temperature,  only  the  White 
cell  data  on  the  pressure-broadened  water  vapor  extinction  coefficient 
are  plotted  in  Figure  21.  These  room-temperature  data  were  for  a  water 
vapor  partial  pressure  of  14.6  torr,  buffered  by  nitrogen  to  a  total 
pressure  of  760  torr.  Also  plotted  for  comparison  in  Figure  21  are 
data  from  Nordstrom  et  al.  [60]  and  McCoy  et  al.  [13],  Peterson  et 
al.  concluded  from  the  strong  quadratic  pressure  of  the  absorption  coef¬ 
ficient  and  the  assumed  Lorentz-line  shape  of  the  far  wings,  that  self¬ 
broadening  (l^O-H^O  interactions)  is  much  more  important  in  its  con¬ 
tribution  to  absorption  at  these  wavelengths  than  is  foreign-broadening 
(H,,0-N2  interactions). 

The  systematic  uncertainties  discussed  in  the  conclusion  to  the 
previous  paragraphs  on  3-5-pm  continuum  are  of  course  present  in  these 
measurements  in  the  8-14-pm  region.  However,  the  characteristic  contin¬ 
uum  water  vapor  absorption  in  this  region  is  roughly  an  order  of  magni¬ 
tude  greater  than  that  in  the  3-5-pm  region.  Thus,  a  conservative  error 
to  associate  with  continuum  data  in  the  8-14-pm  region  is  ±10-15  percent. 

The  proper  model  for  the  temperature  dependency  has  been  a  problem 
for  some  time.  Early  workers  either  neglected  the  temperature  effect 
entirely  or  assumed  the  same  temperature  function  for  both  the  self- 
and  foreign-broadening  terms.  Kunde  and  Maguire  (62]  seem  to  be  the 
first  to  pouit  out  that  the  selr-  and  foreign-broadening  terms  should 
have  different  (indeed  inverse)  temperature  dependencies,  as  shown  in 
Figure  22. 


k(P’PH20'T)  “  kl(T)P  +  k2(T)pH20  (31) 

They  used  the  modified  Van  Vleck-We isskopf  line  shape  for  the  rotational 
liner  et  400  cm  and  the  self-broadened  Lorentz  line  for  the  lines  of 
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Figure  21.  Extinction  coefCicltt.es  from  Peterson  et  al.161], 
Nordstrom  et  al.[60],  and  McCoy  et  al.[13]. 
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Figure  22.  Water  vapor  continuum  absorption  coefficient  for 
the  400-l400-cm_l  region  [62]. 

1595-*cra  H^O  band  at  1400  cm  The  k^  con.ponent  is  normalized  in  a 
similar  fashion  to  the  foreign-broadened  water  vapor  components  at  400 
and  1400  era  These  tenperature  dependencies  were  indicated  by  Tuer 
[63]  as  shown  in  Figure  23,  where  LOWTRAN  II  had  no  temperature  depen¬ 
dence  (except  through  the  density  variation  at  constant  pressure);  and 
Roberts  et  al,  [64  j  suggested  using  the  same  functional  dependence  for 
and  k2» 
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Figure  23.  Comparison  of  various  8-12- ym  water  continuum 
models  [water  content  (14  g/ms)  and  tctal 
pressure  (1  atm)  correspond  to  midlatitude 
summer  model]  [63]. 


Very  recently,  Burch  and  Gryvnak  [65]  presented  some  new  data  on 
H20  continuum  absorption  just  beyond  5  Urn  where  the  CO  laser  operates. 

By  working  in  the  wing  just  below  the  strong  6.3  Um  H90  band,  they  in¬ 
vestigated  the  shape  factor,"  x»  for  nitrogen— broadened  versus  self- 
broadened  lines.  They  concluded  that  both  foreign-  and  self-broadened 
H2°  lines  absorb  more  over  a  large  portion  of  the  wings  than  do  Lorentz- 
shaped  lines,  and  that,  in  the  wings,  self-broadened  lines  absorb  more 
than  nitrogen-broadened  lines. 

Besides  the  above  data  and  the  CO^  and  N2  continua  already  men¬ 
tioned  that  were  measured  by  Burch  and  have  been  incorporated  in  the 
AFGL  Laser  Code,  there  have  been  other  continua  discussed  in  the 
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literature.  (For  example,  Ashcheulov  et  al.  [66]  and  Bolle  [67]  measured 
continuum  from  8  to  25  um;  and  Gates  [68]  aeasured  water  vapor  continuum 
from  0.872  to  2.537  um.)  However,  these  continue  measurements  seem  largely 
unsupported,  and  for  conservatism  in  safety- standards  probably  should  be 
neglected  for  present  purposes. 

Molecular  Scattering — Scattering  of  laser  radiation  depends 
strongly  on  the  size  of  the  particle  relative  to  the  wavelength  of  the 
radiation.  When  the  particle  size  is  much  smaller  than  the  wavelength, 
as  it  is  for  molecules,  the  Rayleigh  theory  gives  a  good  representation 
of  the  scattering  process.  The  volume-scattering  coefficient  for  Ray¬ 
leigh  scattering  can  be  expressed  as: 


(32) 


where 

N,V,n  =  The  number  density,  the  volume,  and  refractive  index 
of  the  molecules 

\  =  Wavelength  of  the  radiation 

n^  =  Refractive  index  of  the  medium  in  which  the  molecules 
are  suspended. 

McClatchey  and  d'Agufi  [42]  propose  a  semiempirical  expression 
based  on  the  equation  for  o  ,  obtained  by  fitting  molecular-scattering 
coefficients  published  by  Penndorf  [69],  This  expression  is  given  in 
terms  of  the  pressure  (atm)  and  temperature  (K)  of  the  atmosphere  and 
the  frequency  of  the  laser  v  (cm  ^) : 

a  »  2.677xlO~17  p/Tv®  (33) 

m 

where  m  is  4.0117  and  0  is  in  units  of  km  It  is  interesting  to  note 

m 

that  the  exponent  m  is  slightly  larger  (i.e.,  approximately  0.3  percent) 
than  Rayleigh  theory. 
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It  appears  that  this  expression,  which  is  contained  in  the  AFGL 
LASER  Code,  should  provide  an  adequate  representation  of  the  extinction 
due  to  molecular  scattering  for  the  present  purposes.  This  extinction 
will  be  negligibly  small  for  most  lasers  of  interest,  as  indicated  in 
Figure  24. 

WAVENUMBER  (cm-1) 

50,000  25,000  16,667  12.500  10,000  8000  6667 


Figure  24.  Model  for  extinction  coefficient  due  to  molecular 
scattering  as  a  function  of  laser  frequency  at 
srp  [42]. 
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Laser  Attenuation  Due  to  Atmospheric  Aerosol  Extinction 

In  this  section,  the  basic  electromagnetic  radiation  scattering 
equations  are  reviewed  and  current  predictive  attenuation  models  are 
evaluated.  For  lower  power  laser  radiation,  molecular  and  aerosol  ab¬ 
sorption  and  scattering  are  the  principal  attenuation.  Currently,  the 
aerosol  attenuation  of  laser  radiation  is  so  poorly  predictable  from 
easily  obtained  atmospheric  parameters  such  as  wind  speed,  relative 
humidity,  and  visibility,  that  clear-air  laser  energy  density  limits 
should  be  employed  for  safety  considerations.. 

Theory- — Scattering  is  defined  as  the  process  by  which  a  particle, 
solid  or  liquid,  continuously  extracts  energy  from  an  incident  electro¬ 
magnetic  wave  and  reradiates  that  energy  into  a  solid  angle  centered  at 
the  particle.  Absorption  can  also  occur  whereby  some  of  the  energy  is 
converted  to  heat  and  reradiated  as  broadband  thermal  infrared  radiation. 
The  combined  effects  of  scattering  and  absorption  are  referred  to  as  ex¬ 
tinction.  Scattering  of  an  incident  beam  of  radiation  requires  that  the 
index  of  refraction  of  the  particle  be  different  from  that  of  the  sur¬ 
rounding  medium.  As  one  might  expect,  the  distribution  pattern  of  scat¬ 
tered  radiation  strongly  depends  on  the  ratio  of  the  particle  size  or 
radius  to  the  wavelength  of  the  incident  radiation  (see  Figure  25). 

The  three  primary  parameters  of  scattering  phenomena  are,  therefore, 
the  refractive  index  and  radius  of  the  particle  and  the  wavelength  of 
the  incident  radiation.  The  scattering  of  radiation  by  particles  much 
smaller  than  the  wavelength  of  the  incident  radiation  (air  molecules, 
for  instance)  is  termed  Rayleigh  scattering,  while  the  scattering  of 
radiation  by  particles  of  radius/wavelength  ratios  of  approximately  0.1 
and  greater  is  referred  to  as  Mie  scattering. 

-A 

Since  atmospheric  particle  sizes  cover  a  range  from  10  pm  to 
10  mm  (Table  16),  and  the  radiation  wavelength  of  interest  is  typically 
from  approximately  0.3  pm  to  20  pm,  both  types  of  scattering  mechanisms 
are,  in  general,  present. 
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S  1 

INCIDENT 

BEAM 

SIZE:  SMALLER  THAN  1/10  THE 
WAVELENGTH  OF  LIGHT 
DESCRIPTION:  SYMMETRIC 

SIZE:  APPROXIMATELY  1/4  THE 
WAVELENGTH  OF  LIGHT 
DESCRIPTION:  SCATTERING  CON¬ 
CENTRATED  IN 
FORWARD  DIRECTION 

SIZE:  LARGER  THAN  THE  WAVELENGTH  OF  LIGHT 

DESCRIPTION:  EXTREME  CONCENTRATION  OF  SCATTERING  IN  FORWARD 
DIRECTION;  DEVELOPMENT  OF  MAXIMA  AND  MINIMA  OF 
SCATTERING  AT  WIDER  ANGLES 


Figure  25.  Angular  patterns  of  scattered  intensity  from 

particles  of  three  sizes:  (a)  small  particles, 
(b)  large  particles,  and  (c)  larger  particles 
[70]. 


TABLE  16.  PARTICLES  RESPONSIBLE  FOR  ATMOSPHERIC 
SCATTERING  [70] 


Type 

Radius  (pm) 

-3 

Concentration  (cm  ) 

Air  molecule 

io-4 

io19 

Aitken  nucleus 

io-3  -  10-2 

io2  -  io4 

Haze  particle'5 

-2 

10-1 

10  -  io3 

Fog  droplet 

1  -  10 

10  -  100 

Cloud  droplet 

1  -  10 

10  -  300 

Raindrop 

io2  -  io4 

-5  -2 

10  -10 

n 

Salt,  dust,  pollen,  etc. 

'’Small  fog  droplet. 
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The  general  approach  to  analytically  modeling  scattering  phenomena 
in  the  atmosphere  has  been  to: 

•  Calculate  the  spectral  scattering  due  to  a  single  particle 
of  radius  r  and  index  of  refraction  n. 

•  Determine  (or  postulate)  a  size  distribution  of  particles 
in  the  atmosphere  as  a  function  of  altitude. 

•  Calculate  the  aggregate  effect  of  all  the  particles  in  a 
particular  radiation  path. 

The  attenuation  of  small-diameter  laser  beams  can,  in  general,  be  deter¬ 
mined  to  sufficient  accuracy  by  considering  only  single  scattering  ef¬ 
fects.  Multiple  scattering  effects  are  much  more  difficult  to  model  and 
calculate. 

The  basic  equation  for  the  scattering  process  begins  with  a  defi¬ 
nition  for  the  angular  scattering  cross  section  of  a  single  particle: 


v(9) 


1A(0)  dX 

-E^dT 


cm2-sr-1 


(34) 


where 


°pX<0> 


angular  scattering  cross  section  for  a  single 
particle  at  wavelength  X 


1(0)  ■  spectral  radiant  intensity  at  wavelength  X 

and  angle  (0)  -  symmetry  about  the  propagation 
axis  is  assumed  (w/sr-1) 

E.  ■  spectral  irradiance  at  wavelength  X  (w/cm-2). 

A 


The  total  scattering  cross  section  for  a  single  particle  is  given  by 
the  integral  of  the  angular  scattering  cross  section  over  4t r  sr: 


0px(9)  du  cn2 


(35) 
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where  dw  =  element  of  solid  angle.  0  ^  can  be  thought  of  as  the  equiv¬ 
alent  area  of  an  elemental  surface  that  completely  scatters  a  portion  of 
the  incident  irradiance  E^dA. 

If  the  scattering  aerosol  consists  of  identical  particles  uniform- 

3 

ly  distributed  with  a  density  of  N  particles/cm  ,  then  a  scattering  co¬ 
efficient  can  be  expressed  by: 

8,  =  No  .  era  1  .  (36) 

A  pA 


8^  can  be  regarded  as  a  coefficient  that  expresses  the  fraction  of  a 
unit  area  that  effectively  scatters  an  incident  irradiance  E^dA  per  unit 
pathlength.  The  primary  assumption  here  is  that  the  dimensions  are 
chosen  such  that  the  medium  is  optically  thin  in  a  unit  pathlength  (i.e., 
6  £  0.2). 


One  might  intuitively  speculate  that  the  maximum  value  o  ,  can 
2 

take  on  is  fir  ,  the  geometric  cross-sectional  area  of  a  particle  of 


radius  r.  Such  is  not  the  case.  An  efficiency  factor  (Q^)  is  defined 
as  tha  ratio  of  the  effective  cross-sectional  area  to  the  geometric 


cross-sectional  area: 


o 

2 

TIT 


(37) 


This  function  is  strongly  dependent  on  the  ratio  of  the  particle  radius 
to  the  wavelength  (A)  of  the  scattered  radiation  (Figure  26).  As  can 
be  seen,  the  function  Q  takes  on  a  maximum  value,  near  4,  for  r/A  ratios 
near  1.  The  feature  that  is  particularly  interesting,  however,  is  the 
rapid  increase  of  this  efficiency  factor  from  r/A  values  of  less  than 
0.5  to  r/A  values  near  1.  This  means  that  the  attenuation  due  to  scat¬ 
tering  can  vary  by  a  factor  of  5  to  10  for  a  relatively  small  increase 
in  particle  diameter.  This  is  particularly  important  under  growing  fog 
conditions. 
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Figure  26.  Scattering  efficiency  factor  versus  size 
parameter  for  water  droplets  [70]. 


For  polydispersions  of  scattering  aerosol,  can  be  expressed  as 
an  integral  over  a  size-distribution  density  function  n(r)  multiplied 
by  the  total  cross-sectional  area: 


n(r)o^(r)  dr 


This  scattering  coefficient  can  also  be  expressed  in  terms  of  the  effi¬ 
ciency  factor  : 


00 

-  f  nr* 

»  A 


Q^(r)n(r)  dr 


In  general,  only  a  range  of  particle  radii  rj-^  interest;  there¬ 

fore: 

r2 

Bx  -  T 1  j  Qx(r)n(r)  dr  (40) 

ri 

where  n(r)  is  the  number  of  particles  per  unit  volume  and  per  unit  inter¬ 
val  of  radius.  The  total  number  of  particles  (N)  in  a  unit  volume  is 
given  by: 


N 


n(r)  dr 


(41) 


Several  particle-size  distribution  models  n(r)  have  been  proposed 
to  describe  observed  distributions  analytically.  The  most  generalized, 
currently  popular,  distribution  model  is  one  proposed  by  Deirmendjian 
[71]  and  is  referred  to  as  the  Deirmendjian  model.  This  analytical 
model  has  the  form: 


,  .  at  -br' 
n(o  »  ar  e 


where  r  is  the  particle  radius  and  a,  b,  a,  and  y  are  positive  constants. 

For  a  particular  choice  of  a  and  y,  the  other  two  constants,  a  and 
b,  can  be  determined  from  a  count  of  the  total  number  of  particles  in  a 
unit  volume,  and  the  mode  radius  rmocje»  where  the  particle  concentration 
is  a  maximum  (i.e.,  curve  derivative  is  zero).  Then: 


b 


tt 


yr 


Y 

mode 


(43) 
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and 


a 


b 


_ _ 

-(a  +  l)  r  a  +  1 

Y  Y 


where  T  is  a  gamma  function.  Typical  haze  distributions  predicted  by 
this  model  are  shown  in  Figure  27. 


Figure  27.  Particle-size  distributions  for  various  haze 

models  [71]:  H  —  high-altitude  stratospheric 
dust  particles;  L  —  continental  aerosols;  and 
M  —  maritime  types  of  aerosols. 
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A  second,  less  general  size  distribution  function  was  proposed 
by  Junge  [72],  and  is  referred  to  as  a  power  law  size  distribution  func¬ 
tion.  This  model  has  the  form: 


n(r)  =*  cr 


dN 

d ( log  r) 


where  c  and  v  are  constants.  The  model  assumes  a  power  law  relationship 
between  the  change  in  particle  concentration  and  the  change  in  the  log 
of  the  radius,  or  is  valid  for  only  the  straight-line  portion  of  the 
[dN/d(log  r)]  versus  log  (r)  curve  (Figure  28). 


LOG  RADIUS  (pm) 


Figure  28.  Particle- size  distribution  showing  range 
of  the  power  law  relationship  [72]. 
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The  constant  (-v)  is  determined  from  the  slope  of  the  straight- 
line  portion  of  the  curve,  and  C  c«n  be  determined  from  a  total  particle 
count  in  a  unit  volume  (N)  by  the  following  equation: 


C 


_ -vN 

0.434(r2"v  - 


(46) 


Current  Prediction  Models — A  number  of  current  prediction  models 
for  aerosol  extinction  were  evaluated  with  regard  to  their  application 
to  low-power  laser  transmission  through  the  atmosphere.  The  salient 
features  of  the  models  are  briefly  reviewed  here.  Published  sources  of 
the  model  descriptions  are  also  included. 

The  overwhelming  evidence  from  the  OPAQUE  measurement  program  is 
that  the  aerosol  extinction  coefficients  cannot  be  predicted  with  ac¬ 
ceptable  accuracy  by  using  easily  measured  atmospheric  parameters  as 
inputs.  Since  all  of  the  models  reviewed  use  the  easily  measured 
parameters,  such  as  relative  humidity,  wind  speed,  and  visibility,  as 
inputs,  none  of  them  would  be  adequate  predictive  models  from  a  safety 
standpoint.  The  six  models  reviewed  are  referred  to  as  Bamhardt  and 
Streete,  Wells-Ga l-Munn ,  Hanel,  LOWTRAN  IV,  Roberts,  and  Science  Appli¬ 
cations. 

The  Barnhardt  and  Streete  model  [73]  was  one  of  the  early  attempts 
to  predict  the  effects  of  aerosol  scattering  on  the  optical  transmittance 
in  the  infrared.  The  aerosol  size  distribution  models  used  were  a  com¬ 
bination  of  a  continental  model  by  .Junge: 


dn(r) 

dr 


0.434  Cr 


-(Y+l) 


(47) 


and  a  maritime  model  by  Deirmendj Ian: 


dn(r)  0,434  L  -ar 

— - -  m  -  be 

dr  r 


(48) 
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where  n(r)  is  the  particle  concentration  in  particle/cm  ,  r  is  the  par¬ 
ticle  radius,  and  a,  b,  c,  and  y  are  the  constants. 

The  relative  humidity  is  included  in  a  growth  factor  (F)  with  the 
following  assumptions: 

•  A  particle  of  any  radius  will  undergo  the  same  fractional 
growth  for  the  same  change  in  relative  humidity. 

•  All  particles  grow  at  the  same  rate. 

•  The  dry  nuclei  have  an  index  of  refraction  of  n=1.54. 

•  The  dilute  water  droplet  is  assumed  to  have  a  refractive 
index  of  n=1.33. 

The  refractive  index  is  assumed  to  vary  between  these  values  as  a  func¬ 
tion  of  the  relative  humidity: 

n  =  1.54  4  0.03  P.n(l-RH)  .  (49) 

The  growth  factor  is  given  by  the  expressions: 

r  -  rQF  (50) 


F 


(51) 


where  r^  is  the  particle  radius  of  0  percent  relative  humidity,  and  3 
is  the  constant  for  a  particular  species  of  condensation  nuclei.  The 
scattering  coefficients  were  determined  by  assuming  a  particular  combi¬ 
nation  of  continental  and  maritime  aerosols  and  relative  humidity  and 
then  calculating,  with  a  digital  computer  program,  the  scattering  coef¬ 
ficients  as  a  function  of  wavelength,  using  Mie  scattering  theory. 

This  Barnhardt  and  Streete  model  does  not  apply  to  atmospheric 
conditions  where  the  relative  humidity  is  approaching  100  percent  or 
for  supersaturated  cases.  It  therefore  has  limited  utility  for  low- 
visibility  or  fog-formation  conditions. 
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The  Wells-Gal-Munn  (WGM)  maritime  aerosol  attenuation  model  [74] 
generates  a  particle-size  distribution  that  is  a  function  of  particle 
radius  (r)f  relative  humidity  (RH),  altitude  (h),  visibility  (v) ,  and 
wind  velocity  (U)  . 

The  overall  particle— size  distribution  model  is  a  combination  of 
the  Deirmendjian  size  distribution  model  for  maritime  aerosols  and 
Junge's  size  distribution  model  for  continental  aerosols. 

The  Deirmendjian  size  distribution  model  used  for  the  maritime 
aerosols  was  of  the  form: 


v 

/  \  H  /  r*)\ 

n(r)  =  ar  e  (52) 

where  r  is  the  particle  radius  for  the  original  Deiiaendjian  model,  and 
a,  b,  i,  and  y  were  constants.  In  the  WCM  model,  a  and  y  are  functions 
of  the  wind  velocity  (U): 


a  =  CA  +  C2l.'d  (53) 

1  =  d  +  pUq  .  (54) 

An  aerosol  growth  factor  (f)  is  included  that  is  a  function  of  relative 
humidity  (RH) : 


F  »  1  -  0.9  Ml-RH)  .  (55) 

The  Junge  size  distribution  model  is  used  for  the  continental 
aerosols : 


n(r) 


(56) 


A  normalization  parameter  is  included  to  define  low-wind  conditions, 
and  a  mixture  parameter  is  included  to  define  the  ratio  of  maritime/ 
continental  aerosols  present  in  the  atmosphere  under  consideration. 
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Exponential  scale  height  multipliers  are  used  to  describe  the 
particle-size  distributions  as  functions  of  altitude.  These  are 
selected  from  a  table  look-up.  Separate  (different)  tables  are  pro¬ 
vided  for  the  continental  and  maritime  distribution. 

Finally,  the  function  dependence  on  visibility  (V)  is  included 
as  multiplier  $  that  is  inversely  proportional  to  V: 


0  = 


3.91 

K0.55V  • 


(57) 


A  Mie  scattering  computer  code  is  then  used  to  determine  the  extinction 
coefficients  as  functions  of  wavelength.  The  index  of  refraction  of 
the  particles  is  based  on  the  assumption  that  all  of  the  particle  nu¬ 
clei  are  NaCl. 

While  this  model  does  consider  wind  velocity  as  a  factor,  the 
fact  that  the  assumed  coefficients  have  a  discontinuity  at  a  particular 
wind  velocity  is  an  indication  that  the  relationship  between  particle 
density  and  wind  velocity  is  complex  and  not  easily  predictable. 

The  decrease  in  particle  density  and  optical  attenuation  with 
altitude  is  not  a  realistic  model  for  conditions  of  low-hanging  clouds 
and  certain  types  of  fogs. 

The  growth  of  atmospheric  models  as  a  function  of  nuclei  size, 
chemistry,  relative  humidity,  and  previous  history  is  an  exceedingly 
complex  phenomenon.  Therefore,  the  simple  relationships  assumed  are 
subject  to  substantial  uncertainties  in  low-visibility,  high-relative- 
humidity  situations.  Further,  as  will  be  demonstrated  later  in  this 
report,  large  variations  in  infrared  extinction  can  occur  with  regard 
to  constant  A=0.55  visible  ranges. 

The  Hanel  Model  175]  provides  a  theoretical  evaluatior  of  aerosol 
particle  mass,  size,  mean  density,  and  mean  refractive  index  as  func¬ 
tions  of  relative  humidity.  The  model  is  rigorous  and  considers  the 
curvature  (radius)  and  dissolved  nucleat ion-site  material  in  the  aero¬ 
sol  particle,  and  the  equilibrium  partial  vapor  pressure  over  the 
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particle.  Particle  growth  (or  shrinkage)  as  a  function  of  particle  radi¬ 
us  and  dissolved  salts  is  considered  and  discussed  in  detail.  This 
riiodel,  however,  has  its  primary  application  to  relatively  stable  light- 
haze  conditions,  since  thermodynamic  equilibrium  was  assumed  in  the 
analytic  derivations.  The  nonequilibrium  conditions  associated  with 
the  formation  or  disintegration  of  clouds  or  unstable  fogs  were  not  con¬ 
sidered.  The  OPAQUE  measurement  data  [76,77]  and  the  laser  transmission 
data  of  Chu  and  Hogg  [78]  have  shown  optical  transmittance  variations  of 
two  to  three  orders  of  magnitude  for  essentially  the  same  humidity  con¬ 
ditions  (RUM).  The  theoretical  aerosol  growth  modeling  of  Neiburger 
and  Chien  [79]  has  also  shown  that  the  particle  concentrations  in  the 
1-1 0-ytn -radius  range  can  increase  by  two  orders  of  magnitude  in  five- 
minute  periods  under  nonequilibrium,  saturation  conditions. 

The  AFGL  LOWTRAN  series  of  atmospheric  transmittance  codes  were 
generated  initially  to  predict  molecular  absorption  effects  on  the  at¬ 
mospheric  optical  transmittance  and  radiance  in  the  spectral  region 
from  0.25  to  28.5  pm.  The  inclusion  of  scattering  effects  on  optical 
transmittance  has  been  a  recent  addition  [80]. 

The  spectral  volume  aerosol  extinction  and  absorption  coefficients 
are  calculated  for  several  "standard"  aerosol  atmospheres  (maritime, 
continental,  urban,  and  rural)  for  visual  ranges  of  23  km  and  5  km. 

The  spectral  aerosol  extinction  and' scattering  coefficients  for  an  ar¬ 
bitrary  visual  range  (V)  are  determined  by  simple  interpolation: 

(I  -  -L) 

o(X,V)  -  Ofx.vo  +  [o(X,V2)  -  o(X,V2)]  (58) 

(v7  "  v7/ 

where  is  the  23-km  visual  range  and  is  the  5-km  visual  range. 

A  functional  dependence  is  also  provided  for  an  altitude  correc¬ 
tion.  This  model  assumes  the  same  spectral  dependence  for  all  aerosol 
distributions  included.  The  model  will  provide  calculations  for  visu¬ 
al  ranges  less  than  5  km;  the  results,  however,  should  be  applied  with 
caution. 
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The  model  apparently  is  not  meant  to  be  used  under  low-visibility 
fog  conditions.  Specific  aerosol  distributions  cannot  be  included  in 
the  model,  or  specific  dependence  on  relative  humidity.  Since  the 
LOWTRAN  TV  model  code  is,  by  far,  the  most  popular  atmospheric  trans¬ 
mission  predictive  code  used  currently,  it  is  unfortunate  that  a  more 
realistic  low-visibility  aerosol  predictive  capability  is  not  included. 

The  R.  E.  Roberts  (Institute  for  Defense  Analysis)  attenuation 
model  [81]  is  based  on  a  simple  phenomenological  approach  to  aerosol 
optical  extinction  modeling.  Basically,  the  model  correlates  the  spec¬ 
tral  extinction  coefficients  due  to  aerosols  with  the  relative  volume 
of  particulates  in  an  atmospheric  path.  Sev^^al  measured  and  theoret¬ 
ical  particle-size  distributions  were  integrated  to  determine  the  rela¬ 
tive  water  content.  The  attenuation  coefficients  at  several  wavelengths 
in  the  visible  and  infrared  were  calculated  for  these  particle-size  dis¬ 
tributions,  using  Mie  scattering  theory.  The  calculated  extinction  co¬ 
efficients  were  then  related  to  the  relative  water  content.  Log-log 

plots  of  the  extinction  coefficient  (for  A*s  =  1.0,  4.0,  and  10.0  pm) 

-1  3 

in  km  versus  relative  water  content  in  g/m  were  relatively  straight 

lines  and  with  approximate  data  spread  of  less  than  a  factor  of  two. 

For  wavelengths  greater  than  the  particle  radius,  Roberts  sug¬ 
gested  that  the  extinction  cross  section  should  be  proportional  to  the 
radius  cubed,  or 


ext 


"v  r 


(59) 


and  for  particle  radii  much  greater  than  the  wavelength,  Mie  scattering 
would  predict  a  relationship: 


“ext  %  ‘  (60) 

Therefore,  it  was  argued  that  the  relationship  between  the  A  -  .l.O-Dm 
extinction  coefficient  and  the  A  =*  10.0-pm  extinction  coefficient  would 
be  of  the  form: 


6(10.d  pm) 


Ca3 


3/2 

(1.0  pm) 


(61) 
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or  that  the  B.,,.  .  would  Increase  as  a  function  of  relative  water 

(10.0  pm) 

content  at  a  faster  rate  than  the  B^j  ^  ^  extinction  coefficient. 

The  calculated  data  generated  by  Roberts,  in  general,  support  this 
hypothesis,  with  the  plotted  data  showing  a  relationship: 


e  =  c  r1-49 

^(lO.O  pm)  Lr(1.0  pm) 


(62) 


While  this  approach  of  attempting  to  relate  the  optical  attenua¬ 
tion  due  to  aerosols  to  a  parameter  such  as  the  total  water/unit  area 
in  a  given  path  is  an  improvement  over  attempts  to  correlate  attenuation 
with  relative  humidity,  some  problems  still  remain.  Specifically: 

The  problem  of  determining  the  water/particulate  content 
in  a  path.  Direct  measurements  are  not  easy. 

*  If  visibility  is  used  as  a  parameter  to  determine  water 
content,  then  the  model  is  subject  to  all  previously 
discussed  pitfalls  with  regard  to  using  visible  radia¬ 
tion  attenuation  as  a  measure  of  infrared  attenuation. 

'  The  OPAQUE  measurement  data  [  / 7  J  indicate  that  large 
variations  can  be  expected  in  attempting  to  correlate 
visible  extinction  coefficients  (X  =  0.55  pm)  with  the 
infrared  window  band  extinction  coefficients  (AX  = 

3-5  and  8-12  pm). 

The  SAI  aerosol  attenuation  model  [82]  uses  two  basic  types  of 
input  data:  (1)  the  ground-level  visibility,  and  (2)  the  aerosol-mass 
content.  This  approach  has  the  unique  feature  of  not  being  functionally 
dependent  on  relative  humidity.  The  relative  humidity,  while  rather 
easily  measured  at  lower  values,  becomes  increasingly  difficult  to  ac¬ 
curately  measure  as  it  approaches  one  and  reaches  supersaturation  values. 
Further,  it  has  been  shown,  theoretically  and  experimentally,  that  rela¬ 
tive  humidity  is  a  poor  indicator  of  infrared  extinction  under  low-vis¬ 
ibility  conditions.  However,  although  the  aerosol-mass  content  may  be 
a  more  useful  parameter  for  predicting  infrared  extinction  accurately, 
determining  the  mass  content  from  measurable  parameters  is  not  simple 
and  straightforward. 

This  aerosol-mass-content  parameter  is  determined  in  the  SAI  model 
by  calculating  the  Mie  scattering  extinction  coefficients  as  a  function 
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of  wavelength  for  measured  particle-size  distribution.  The  size  dis¬ 
tributions  are  integrated  to  determine  the  aerosol-mass  content  (as¬ 
sumed  to  be  water).  The  volume  optical  extinction  coefficients  are 
log-log  plotted  as  a  function  of  the  aerosol  mass  (liquid  water  con- 
tent)  in  g/m  .  From  these  plots,  the  slopes  and  intercepts  are  deter¬ 
mined  for  analytically  correlating  the  extinction  coefficients  with 
liquid  water  content.  The  visible  volume  extinction  coefficient  (A=0.55 
lim)  is  related  to  the  surface  horizontal  visibility  range  by  the  inverse 
relationship: 


(1*0.55  um) 


=  3.912  V 


(63) 


Therefore,  the  visibility  (V)  is  used  to  determine  the  surface  aerosol- 
mass/water  content. 

Measurements  were  made  at  Grafenwohr,  West  Germany,  to  determine 
the  vertical  profile  of  the  aerosol  mass.  In  the  SAI  model,  the  verti¬ 
cal  profile  is  approximated  by  exponential  functions  of  the  form 


(z) 


W  e 

"o 


pz 


(64) 


where  z  is  the  vertical  height  above  Earth's  surface.  Changes  in  che 
profile  due  to  increasing  or  decreasing  aerosol-mass  density  are  han¬ 
dled  by  discrete  changes  in  the  factor  p.  The  total  liquid  water  per 
unit  along  a  slant  path  is  determined  by  integrating  the  altitude  den¬ 
sity  profile  along  the  slant  path.  From  this  information,  the  aerosol 
extinction  coefficients  are  determined  and  the  spectral  transmittance 
functions  are  calculated. 

This  SAI  aerosol  attenuation  model  is  an  improvement  over  previ¬ 
ous  models  that  used  relative  humidity  as  an  input  parameter,  and  is 
best  suited  for  predicting  attenuations  under  hazy  conditions.  However, 
this  model  has  substantial  deficiencies  with  regard  to  accurate  predic¬ 
tion  of  infrared  attenuation  under  low-visibility  conditions.  Determi¬ 
nation  of  the  infrared  attenuation  is  functionally  dependent  on  the 
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visible  range  (at  A=0.55  Pm)  through  the  aerosol-mass  content.  OPAQUE 
data  [77]  and  the  laser  transmittance  data  of  Chu  and  Hogg  [ 78 1  have 
shown  that  the  visibility  can  vary  by  two  orders  of  magnitude  for  a 
relatively  constant  IR  transmittance.  Further,  in  the  model  the  sur¬ 
face  visibility  and  the  vertical  aerosol— mass  density  in  mass/unit  area 
(referred  to  as  the  liquid  water  column  density)  are  regarded  as  inde¬ 
pendent  input  variables.  Smith's  comment  [82 1  that  "there  is  prob¬ 

ably  some  correlation  between  the  two  in  actual  situations"  appears  to 
be  a  gross  understatement. 

Measurement  Results — A  number  of  predictive  models  for  infrared 
attenuation  have  attempted  to  correlate  radiative  extinction  in  the  two 
primary  atmospheric  window  regions  (3-5  pm  and  8-14  pm)  with  visibility 
(i.e.,  transmission  in  the  visible  portion  of  the  spectrum).  The  haz¬ 
ards  in  attempting  such  correlations  were  demonstrated  by  some  very 
careful  measurements  by  Chu  and  iiogg  [78]  (Figure  29).  In  the  measure¬ 
ment  data  shown,  the  visibility  (0.63  pm)  radiation  degenerates  drasti¬ 
cally  (more  than  three  orders  of  magnitude)  as  the  fog  becomes  denser, 
while  the  10.6-pm  transmittance  is  relatively  unattenuated  and  the  3.5-- 
pm  radiation  is  only  slightly  decreased.  The  upper  curves  show  the 
transmission  in  the  3.5-and  0.63-pm  regions  improving  as  a  fog  dissi¬ 
pates.  As  might  be  expected,  as  the  fog  droplets  shrink  in  size,  the 
3-5 -pm  transmission  improves  first  as  the  mean  droplet  radius  passes 
through  3.5  pm.  The  visible  transmittance  then  improves  drastically 
after  the  mean  radius  shrinks  below  0.6  pm. 

Theoretical  work  by  Neiburger  and  Chien  [79]  has  predicted  rapid 
changes  in  the  size  distribution  function  for  increasing  fog  conditions. 
A  change  of  two  orders  of  magnitude  in  the  concentration  at  5  pm  over  a 
five-minute  interval  is  shown  (Figure  30).  A  biaodal  distribution  is 
also  predicted  for  certain  development  stages.  Since  fog  development 
occurs  under  essentially  saturated  or  supersaturated  conditions  (RH  « 
100  percent),  large  changes  in  the  visible  and  near- infrared  transmit¬ 
tance  can  occur  for  essentially  constant  high-relative-humidity  condi¬ 
tions.  This  variability  is  graphically  demonstrated  by  the  results  of 
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Figure  29.  Measurement  of  2.6-km  transmission  loss 
in  light  fog;  O-dB-signal  level  in  clear 
weather  [78]. 

the  extensive  OPAQUE  measurement  program  [77]  (Figure  31).  As  shown, 
the  visible  extinction  coefficient  (0.55  Urn)  varies  by  nearly  two  orders 
of  magnitude  for  high  humidity  conditions.  The  difficulties  in  attempt¬ 
ing  to  use  visibility  in  predicting  3-5-pm  and  8-12-pm  extinction  coef¬ 
ficients  are  also  demonstrated  by  the  OPAQUE  data  (Figure  32).  As  the 
data  show,  it  is  not  possible  to  reliably  predict  the  infrared  extinc¬ 
tion  coefficient  to  within  a  one-half  order  of  magnitude  given  visible 
extinction  values. 
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PARTICLE  DROPLET  RADIUS  (pm) 

Figure  30.  Theoretical  size  distribution  of  fog  droplets 
as  a  function  of  elapsed  time.  The  time  in 
minutes  after  start  is  shown  for  each  curve 
f7y  1. 
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Figure  31.  Variation  of  visible-aerosol  extinction  coef¬ 
ficient  with  relative  humidity.  Derived  from 
OPAQUE  Netherlands  3/77-5/77  data  [77], 
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AEROSOL  EXTINCTION  COEFFICIENT  AT  »-6  fim  (km' 


Figure  32.  Comparison  of  various  aerosol  models  with  OPAQUE 
dara.  Aerosol  extinction  coefficient  at  0.55  Urn 
versus  aerosol  extinction  coefficient  at  3-5  tun. 
Derived  from  OPAQUE  Netherlands  3/77-5/77  data 
l  771. 
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Atmospheric  Turbulence  Effects 

Turbulence  in  the  atmosphere  can  cause  random  variations  in  the 
optical  properties  of  the  atmosphere,  which  affect  the  statistical  atmo¬ 
spheric  transmission  characteristics  for  laser  radiation.  This  effect 
has  long  been  a  limiting  factor  in  the  resolution  achieved  by  ground- 
based  astronomical  observatories,  and  is  commonly  referred  to  as  atmo¬ 
spheric  scintillation  or  twinkling.  Laser  radiation  can  be  similarly  in¬ 
fluenced  by  atmospheric  turbulence  in  the  path  between  the  transmitter 
and  the  receiver,  with  random  temporal  and  spatial  variations  in  the  re- 
ceived  irradiance.  With  regard  to  safety  considerations,  these  statisti¬ 
cal  variations  imply  the  possibility  of  local  bright  spots  whose  irradi¬ 
ance  can  be  as  much  as  five  times  larger  than  the  average.  In  addition, 
the  turbulence  can  cause  random  spreading  of  the  laser  beam  that  would 
reduce  the  irradiance  (per  unit  area  of  collector).  For  these  reasons, 
it  is  important  to  investigate  both  of  these  atmospheric  turbulence  ef- 
*  fects  on  the  transmission  of  laser  radiation,  and  to  determine  the  best 
models  for  predicting  the  devel  of  atmospheric  turbulence  and  its 
effects  on  laser  transmission. 

The  level  of  turbulence  in  the  atmosphere  typically  varies  strongly 
with  time  of  day,  wind,  cloud  cover,  local  terrain  features,  and  altitude. 
The  discussion  in  the  following  paragraphs  evaluates  the  available  infor¬ 
mation  on  the  theories,  experiments,  and  models  pertinent  to  the  expected 
turbulence  levels  in  the  atmosphere  as  a  function  of  these  and  other 
pertinent  parameters.  The  section  entitled  "Turbulence  Effects"  presents 
an  evaluation  of  the  available  information  on  the  effects  of  such  turbu¬ 
lence  that  are  important  to  personnel  safety.  This  evaluation  also 
addresses  the  theories,  experiments,  and  models  relevant  to  these  effects. 

Atmospheric  Turbulence  Levels — Turbulence  in  the  atmosphere  is  the 
result  of  several  mechanisms  such  as  wind  shear  over  terrain,  vegetation, 
and  buildings,  and  the  result  of  thermal  nonequilibrium  between  various 
layers  in  the  atmosphere  or  with  the  ground.  Transient  and/or  spatially 
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norumifonn  heating  of  the  ground  due  to  the  sun  (i.e.,  sunrise,  sunset, 
and  broken  cloud  cover)  is  the  main  source  of  the  local  nonequilibrium 
between  the  ground  and  the  air  just  above  it. 

In  these  paragraphs,  the  theory,  measurements,  and  models  overlap 
so  much  that  they  were  not  called  out  in  separate  headings.  However,  as 
much  as  possible,  the  topic  is  addressed  in  that  order:  a  general  histor¬ 
ical  discussion  of  turbulence  comes  first,  followed  by  specific  theory  and 
measurements  on  atmospheric  turbulence.  Low-altitude  turbulence  will  be 
discussed  first  and  then  high-altitude  turbulence. 

Turbulence  in  a  fluid  has  long  been  the  subject  of  theoretical  and 
experimental  investigations.  The  origin  of  turbulence  and  the  mechanism 
affecting  the  transition  from  laminar  to  turbulent  flow  are  of  fundamental 
importance  for  the  entire  field  of  fluid  mechanics.  Fluid  dynamists  tra¬ 
ditionally  characterize  the  flow  of  gases  by  the  Reynolds  number,  which 
can  be  thought  of  as  the  ratio  of  the  inertial  forces  to  viscous  forces 
acting  on  the  fluid.  They  have  found,  in  controlled  experiments,  chat 
there  is  a  sudden  transition  from  laminar  to  turbulent  conditions  as  the 
Reynolds  number  passes  some  critical  number  (approximately  2300). 

In  the  open  atmosphere  over  the  ground  there  are  a  number  of  com¬ 
plicating  factors,  such  as  the  ground  heating  or  cooling  faster  than  the 
air,  which  produce  convective  currents.  These  unstable  conditions  (some¬ 
times  called  active  conditions)  are  frequently  quantified  by  the 
Richardson  number: 
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The  gradient  in  0  is  given  as: 
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(66) 


_3 

where  y  is  the  adiabatic  iapse  rate  of  a  parcel  of  dry  air  (9.8x10  K/m). 

Unstable  conditions  are  characterized  by  negative  values  of  R^,  and  stable 
conditions  by  positive  values.  Nighttime  conditions  are  usually  stable, 
since  convection  is  inhibited  by  warmer  air  above  the  cooler  surface  air. 

The  parameter  generally  used  as  a  measure  of  the  level  of  turbu- 

2 

lence  is  the  temperature  structure  constant  (C^,),  which  is  related  to  the 
expected  value  of  the  square  of  the  temperature  difference  between  two 
points  (r^  and  r^)  separated  by  a  small  distance  (r): 

cr  :  r~2/3{[T<r2)  '  l'(r1)]2}  *  (67) 

2 

Wyngaard  and  Izumi  £83]  developed  an  expression  for  C“  near  the  ground  in 
terms  of  a  function  of  : 


C2  =  Z4/3(30/i2)2f(R.) 


(68) 


2 

Measurements  of  C^,  over  a  variety  of  conditions  at  heights  of  5.66,  11.3, 
and  22.6  m,  were  correlated  to  the  Richardson  number,  as  shown  in  Figure 
33.  Note  that  the  turbulence  is  shown  to  become  very  small  for  R^  >  0.2. 
Turbulence  does  develop  even  under  these  highly  stratified  conditions, 
although  at  present  there  is  no  adequate  treatment  for  these  situations. 

The  variation  of  cf  with  altitude  was  also  investigated  theoreti¬ 
cally  and  experimentally  by  Wyngaard  and  Izumi  [83].  They  show  semiem- 
pirical  relations  for  CT  in  the  surface  layer  (i.e.,  up  to  approximately 
25  m)  in  terms  of  the  Monin-Obukhov  length  L  and  a  parameter  T*  (see 
Figure  34) .  These  parameters  are  defined  by: 


L  ”U*T/kgQ 


(69) 
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ignre  33. 


HI 

Ihe  dimensionless  temperature-structure 
parameter  versus  Kicliardson  number  [83] 


and 


T*  r  -Q/u*  »  (70) 

where  is  the  kinematic  surface  stress  (i.e.,  the  stress  per  unit  air 

2  2 

density  exerted  by  the  wind  on  the  surface  in  cm  /sec  ),  Q  is  the  surface 
temperature  flux  (in  K  cm/sec),  and  k  is  von  Karman's  constant  (tradition¬ 
ally  taken  to  be  0.4).  The  functional  form  in  this  figure  is  given  by: 

2  -2/3 

CT  =  T*Z  83(Z/L)  (71) 

where 

g3  =  4.9(1  -  7Z/L)~2/3  for  Z/L  <  0 


and 


g3  =  4.9(1  +  2.74Z/L)  for  Z/L  >  0  . 


This  indicates  a  transition  from  a  Z  2  ^  decay  at  low  altitudes  to  a 
-4/3 

Z  decay  at  higher  altitudes  occurring  a  few  meters  above  the  ground. 
-4/3 

The  Z  relationship  has  been  experimentally  substantiated  up  to  alti¬ 
tudes  in  the  100-  to  500-m  range  [84],  Unfortunately,  this  expression 
involves  the  parameters  T^  and  L  which  are  not  easily  obtainable. 
Hufnagel  [85]  recasts  the  expression  for  the  thermally  active  region 
(i.e.,  Z/L  <  0)  into  a  more  convenient  form: 


2xl0-3  Q/+/3 


3  \"2/3 
14000  U*  \ 


hQ 


(72) 


and  gives  approximate  relationships: 

*  0.35z(||)  m/s 


(73) 
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and 


1 


2 

Q  =  Q  sin  T,  -  50  w/m  (74) 

where  s  is  the  solar  zenith  angle  and  is  as  given  in  Table  17  for 
various  conditions. 

TABLE  17.  REPRESENTATIVE  VALUES  OF  Q()  VERSUS 
TERRAIN  AND  CLOUD  COVER  [S5] 


Type  of  terrain 

% 

Coat  L.uousl y 
clear  sky 

% 

Overcast  sky 

Dry  sand  or  lava 

500 

200 

Dry  field  or  brush 

400 

150 

Wet  fields 

200 

70 

There  are  a  large  number  of  measurements  of  the  turbulence  level 
near  the  ground  (see  Table  18).  Several  of  these  measurements  are  rep¬ 
resented  in  Figure  35  in  the  form  of  an  index  of  refraction  structure 
•  2 

"constant,"  C  ,  as  a  function  of  altitude.  This  structure  function  is 
n 

directly  related  to  that  of  temperature  and  is  the  parameter  used  to 
relate  the  turbulence  levels  to  their  effects  on  laser  beam  propagation. 
Considerable  scatter  is  seen  in  the  measured  data  as  might  be  expected 
due  to  the  random  nature  of  the  phenomena.  However,  the  general  levels 
and  trends  of  most  of  the  data  seem  roughly  self-consistent,  with  the 
possible  exception  of  Subramanian 's  measurements  that  seem  high,  and 
Wright  and  Schutz's  minimum  that  seems  low. 

For  comparison  with  models,  the  above  data  were  replotted  in 
Figure  36,  and  the  turbulence  levels  predicted  by  Hufnagel's  model  [85] 
for  extremes  in  solar  heating  and  wind  speed  were  superimposed.  The 
maximum  solar- heat ing  case  (i.e.,  clear  sky,  sun  at  zenith  over  dry  sand 
or  lava)  produces  h igli-turbulence  levels  as  reported  by  Subramanian, 
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TABLE  18.  LOW-ALTITUDE  MEASUREMENTS  OF  REFRACTIVE  INDEX  STRUCTURE  PARAMETER 


10°  101  102  wr*  10* 

HEIGHT  ABOVE  GROUND  LEVEL  (ml 


Figure  35.  Comparison  of  turbulence  models  with 
measurements  -  low  altitude. 
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Figure  36.  Comparison  of  turbulence  models  with 
measurements  -  low  altitude. 


Iml  has  a  much  taster  decay  rate  with  altitude  than  those  measurements. 
Tiie  turbulence  decay  rates  for  both  maximum  and  minimum  illumination  con¬ 
ditions  with  no  wind,  however,  appear  in  general  agreement  with  the  rest 
of  the  measurements  shown.  Wind  reduces  the  turbulence  levels  much  more 
for  ttie  minimum  heating  condition  than  the  maximum,  according  to  this 
model . 

It  should  be  realized  that  the  structure  constant  io  very  sensi¬ 
tive  to  local  winds  and  local  terrain  features,  so  the  large  spread  in 
the  measurements  should  not  be  unexpected.  For  these  same  reasons,  a 
simple  model  such  as  Hufnagel's  cannot  tie  expected  to  give  more  than 
general  indications  of  the  levels  of  turbulence  that  can  be  expected. 

For  safety  considerations,  the  upper  models  (i.e.,  maximum  insolation) 
should  probably  be  used  in  all  situations  to  be  conservative,  since  the 
phenomena  are  not  completely  understood  at  the  present  time  and  detailed 
empirical  verification  of  the  models  is  lacking. 

At  higher  altitudes,  the  turbulence  theory  is  poorly  understood 
and  strongly  influenced  by  the  nonstationarity  of  the  entire  boundary 
layer  [85],  However,  there  have  been  numerous  measurements  (see  Table 
19)  and  several  basically  empirical  models  fit  to  the  data.  Some  of 
these  measurements  and  models  are  compared  in  Figure  37.  The  solid 
curves  on  the  figure  are  various  models  by  Hufnagel  and  by  Fried.  There 
is  a  considerable  spread  in  the  measured  turbulence  levels  at  high  al¬ 
titudes,  but  most  of  the  spread  originates  in  tha  extensive  aircraft 
measurements  by  Norris  [93].  Morris'  data  also  indicate  turbulence 
levels  that  are  approximately  an  order  of  magnitude  larger  than  any  of 
the  other  measurements  or  models,  Such  aircraft-borne  measurements 
would  seem  very  difficult  to  make  and  particularly  sensitive  to  bias. 

A  number  of  models  have  been  developed  for  the  high-altitude  re¬ 
gion.  Hufnagel  and  Stanley  [95]  made  one  of  the  first  attempts  to  semi- 
empirically  model  the  observed  variation  of  the  structure  constant  with 
altitude.  The  expression  they  give  is: 


a 


u2  F.2/3  yV 


(75) 
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TABLE  19.  HIGH-ALTITUDE  MEASUREMENTS  AND  MODELS  OF  REFRACTIVE  INDEX  STRUCTURE  PARAMETER 


Model  includes  wind  end  a 
random  variable  lor  lino 
structure 


MEASUREMENTS: 


□ 


Figure  37.  Comparison  of  turbulence  models 

with  measurements  -  hir>h  altitude. 


where 


a  =  Constant 

E  =  Rat*,  of  energy  per  uuu  mass  dissipated  by  viscous  friction 
y  =  Vertical  gradient  of  the  potential  temperature 
3  =  Average  shear  rate  of  the  wind. 

In  their  paper,  Hufnagel  and  Stanley  {95]  present  average  relationships 

for  tne  parameters  of  Eq.  75  as  a  function  of  altitude  from  various 

sources  that  were  used  in  their  model.  This  model  is  indicated  by  the 

broken  curve  in  Figure  38.  The  solid  curve  was  not  discussed  by  them, 

-1/3 

but  presumably  it  is  a  fit  to  data  of  the  form  Z  as  predicted  at  low 
altitudes  by  Tatarski  {91].  Measured  data  (with  and  without  variability 
bars)  are  from  Tatarski  {91],  Zwang  f 96 ] ,  and  Gossard  {97].  Fried  [98] 
later  provided  the  expression: 

C~  =  h  ^  exp(~h/hQ)  (76) 

with  C^=4.2xl0  8=1/3,  hg=3200  m,  and  given  in  units  of  m 

which  is  intended  to  represent  Hufnagel's  model  [95], 

In  1966,  Hufnagel  updated  this  model  [99]  and  included  the  ability 
to  model  distributed  layers  in  the  atmosphere  (see  Figure  39).  Brookner 
[100]  extended  Eq.  76  to  obtain  an  analytical  expression  of  this  Hufnagel 
model  [99]  (including  the  disturbed  layer): 

h  ^  exp(-h/h  )  +  C+  6(h-hD)  (77) 

n  n^  On  ^ 

+  -13  1/3 

where  is  a  random  quantity  with  a  mean  value  of  4.3x10  m  and  a 

standard  deviation  factor  of  (Zn  2)  greater  than  that;  6  is  the  delta 

function.  This  expression  does  not  adequately  represent  Hufnagel's 

model  over  a  very  large  altitude  regime,  however,  and  so  does  not  seem 

to  be  of  much  value  for  the  present  problem. 
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Figure  38.  Early  mode]  of  temperature  structure  constant 

versus  altitude  by  Hufnage!  and  Stanley  compared 
with  measurements  [95], 
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Figure  39.  Hufnagel's  1966  model  for  index  of  refraction 
structure  constant  C  [99], 
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Hufnagel  later  proposed  another  n;jdel  [101],  valid  between  2'J0  and 
2400  m  above  ground,  which  is  similar  ir  form  to  Brookner's  but  more 
representative  of  nature: 


)  In  ? 

(T  =  2. 2xi0  h  (U/27)  e*p(-h/1000) 

+  10~L6  exp (li/l 500)  er(h,t)  .  (78) 

This  expression  contains  a  factor,  l>,  which  is  .-he  root- mean-squared 
(RMS)  wind  speed  from  5  to  20  km,  and  alro  allows  for  the  generation  of 
fine  structure  by  the  Gaussian  random  variable  r,  with  zero  mean  and  a 
variance  of  / 2 ,  For  those  not  interested  in  the  fine  structure,  the 
random  exponential  factor  can  be  replaced  by  the  constant  e  (i.e., 
2.7183).  Hufnagel  also  provides  a  typical  value  for  U  (i.e.,  RMS  wind 
speed  U=18  m/sec),  lie  presents  a  sample  of  random-generated  atmospheric 
turbulence  function  which  aesthetically,  at  least,  looks  similar  to  real 
data  (see  Figure  40).  Hufnagel  i f» 5 ]  recently  gave  an  equation  equiva¬ 
lent  to  Eq.  78  except  without  tee  random  factor  (which  he  said  could  be 
added  if  desired).  In  this  last  discussion,  however,  he  also  suggests 
that  U  can  be  represented  by  a  Gaussian  random  variable  that  varies  from 
day  to  day,  with  a  mean  of  27  m/sec  and  a  standard  deviation  of  9  m/sec. 

These  more  prominent  models  of  high-altitude  atmospheric  turbulence 
discussed  above  were  compared  in  Figure  37  with  the  previousl"  presented 
measurements.  Hufnagel1 s  1974  high-altitude  turbulence  model  [101]  ap¬ 
pears  to  match  the  measured  data  shown,  except  for  Morris’,  quite  well. 
This  model  can  he  evaluated  in  three  different  ways.  The  first  and 
easiest  is  to  simply  use  a  nominal  average  value  for  0  in  the  5-20-km 
region  (e.g.,  27  m/sec)  (as  was  done  for  the  curve  shown  in  Figure  4oa). 
The  second  method  is  to  use  radiosonde  measurements  of  the  wind  speed 
for  the  actual  conditions  and  altitude  of  Interest.  The  third  produces 
a  random  atmospheric  value  by  using  a  random  sample  from  a  Gaussian  dis¬ 
tribution  of  speeds  with  a  nominal  mean  and  standard  deviation  to  get 
the  day-to-day  variability  (e.g.,  27  and  9  m/sec,  respectively ) .  Such 
random  fluctuations  are  frequently  observed  in  !. igh-resolut ion  measure¬ 
ments  at  altitudes  above  the  first  inversion  (see  Figure  40b).  However, 
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Figure  40a.  Sample  of  atmospheric  turbulence  calculated 
with  Hufnapel's  random  model  (with  V  *  18 
m/sec  )  [101]. 
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Figure  40b.  Measured  atmospheric  turbulence  and  air 
temperature  as  a  function  of  altitude 
over  the  ocean  (85). 


105 


the  third  method  would  no*-  be  the  best  from  a  conservative  safety  view¬ 
point.  The  second  method  is,  of  course,  the  best  of  the  three,  but  the 
first  should  be  satisfactory  if  a  safety  margin  were  added.  This  could 
oe  accomplished  sin, ply  by  multiplying  the  coefficient  of  the  second  term 
(i.e.,  2.7x10  ^)  by  a  factor  (e.g.,  1.1). 

Turbulence  L*~f acts' — Atmospheric  turbulence  in  the  path  of  a  laser 
beam  causes  tenrjr.il  and  spatial  variations  in  its  radiation  field.  The 
variations  prodi-'e  localized  regions  of  intensification  that  can  possi¬ 
bly  enhance  sign . f icantiy  the  probability  of  laser  eye  damage.  Existing 
theories,  measurements,  and  models  relevant  to  these  effects  are  re¬ 
viewed  here  to  assess  the  state  of  information  in  this  field  and  make 
recommendations  for  predictive  models  to  aid  in  setting  comprehensive 
safety  standards. 

As  a  field  of  study,  the  effect  of  turbulence  on  optical  transmis¬ 
sion  has  had  a  relatively  long  and  active  history  due  in  part  to  its  im¬ 
portance  in  astronomical  observations.  In  1952,  Chandrasekhar  [102], 
estimated  the  amplitude  scintillation  caused  by  atmospheric  turbulence, 
using  a  geometric  optics  approach.  In  1955,  Muchmore  and  Wheelan  [103] 
extended  Chandrasekhar's  work  by  calculating  the  variance  and  spectral 
correlation  functions  for  the  amplitude,  phase,  and  angle  of  arrival  of 
the  laser  radiation,  using  the  geometric  optics  formulation.  The  Born 
approximation  (i.e.,  method  of  small  perturbations)  was  used  by  Wheelan 
[104]  to  solve  the  scalar-wave  equation  including  diffraction  effects. 

Obukhov  [105]  in  1953  apparently  first  applied  a  method  proposed  by 
Rytov  [106]  to  consider  atmospheric  turhuience  effects  in  sound  and 
light  propagation;  however,  this  work  did  not  receive  wide  circulation 
until  it  was  translated  in  1960  [ 107,91 j.  The  Rytov  method  Is  similar 
to  that  of  Born  in  that  they  are  both  perturbation  techniques;  however, 

Rytov's  method  includes  multiple  scattering,  while  Born's  method  treats 
only  single  scattering  and  is  more  widely  used  to  correlate  measured 
turbulence  effects.  Rytov's  method  is  much  more  severely  limited  In  the 
range  of  applicability  than  originally  thought. 

i  i 
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Rvtov's  theory  results  in  the  following  equation  relating  the 
variance  in  the  log  of  the  amplitude  of  the  laser  irradiance  (cT-p)  at  a 
point  in  space  to  the  atmospheric  turbulence  structure  function  «&: 


J2  =  0.56(2«/A)7/6 


( 


C2(x)(x/L)5^b 
n  / 


(L-x)^°  dx 


(79) 


wiiere  the  integral  is  carried  out  over  the  entire  path  (i.e.,  x  running 
from  zero  to  L) .  If  the  turbulence  is  uniform  over  the  path,  this  can 
be  integrated  to  yield  (1081: 

->2  =  0. 124  (2i./ \ )  7/6  L11/6  C2  .  (80) 

r  n 

This  relat ionsh ip  ties  been  verified  by  several  experiments  (e.g.,  see 
Figure  41),  within  the  limit  of  the  scatter  in  the  measurements.  The 
set  of  data  by  Johnson  et  al,  [109 j  (Figure  41)  indicates  excellent 
agreement  with  linear  theory  up  to  and  slightly  beyond  log  intensity 
standard  deviations  of  unity.  Their  measurements  were  taken  over  a 
460-m  range,  and  the  refractive  index  structure  constant  was  derived 
from  differential  thermometer  data.  A  data  set  by  Kerr  (110]  (Figure 

4°),  which  is  a  composite  of  measurements  at  different  wavelengths,  in- 

2 

dicates  similar  agreement  with  Rytov's  theory  (°j)  but  seems  to  drop 
away  from  the  linear  behavior  early  (i.e.,  at  lower  turbulence  levels). 

Gracheva  and  Gurvich  in  1965  fill]  were  the  first  to  demonstrate 
experimentally  that  the  variance  In  the  irradiance  does  not  continue  to 
increase  1 Inearly  with  the  level  of  turbulence  as  characterized  by  Cn, 
the  so-called  saturation  effect  (e.g,,  see  Figure  43).  This  was  soon 
verified  by  other  experiments  (12-1 5}.  In  1970,  both  Gracheva  et  al. 

( 1 1 3 J  and  Kerr  [110]  measured  the  supersaturat ion  effect,  where  the 
variance  of  the  irradiance  actually  grows  smaller  as  the  turbulence 
level  is  increased  still  more  beyond  che  linear  region  (see  Figure  44). 
Reviews  of  several  reported  measurements  Indicate  that  the  maximum 
standard  deviation  always  seems  to  lie  below  a  value  of  1.6,  as  shown 
In  Figure  45. 
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Figure  41.  Observed  irradiance  statistic  versus 
measured  structure  constant  (linear 
scintillation  region)  (109]. 
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Figure  42.  Experimental  versus  theoretical  log-amplitude 
variance  for  4880  A,  1.15  ym,  and  10.6  ym 
combined  [110]. 
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Figure  43.  Measurements  of  variance  in  log  intensity  as  a  function  of 
the  turbulence  level,  showing  the  saturation  effect  [112]. 


There  has  been  some  controversy  over  the  type  of  statistical  dis¬ 
tribution  function  that  best  represents  these  irradiance  variations. 

It  is  generally  taken  to  be  log-normally  distributed,  but  some  authors 
have  preferred  the  Rice-Nakagami  distribution.  It  currently  appears 
that  the  consensus  is  with  the  log-normal  distribution.  This  distribu¬ 
tion  appears  to  apply  even  outside  the  range  of  linear  effects  into  the 
saturation  region  [112]. 

A  number  of  empirical  models  have  been  developed  to  predict  the 
functional  behavior  of  the  standard  deviation,  but  there  is  no  consensus 
as  to  the  best.  It  appears  that  the  model  by  Johnson  et  al.[109]  fits 
his  data  very  well  (see  figure  44),  and  seems  to  have  the  general  shape 
of  the  other  measurements  (see  Figure  45): 
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Figure  44.  Observed  relationship  between  the  theoretical  stan  - 
dard  deviation  (using  Linear  theory)  and  measured 
standard  deviation;  shows  the  effect  of  supersatura¬ 
tion  [109]. 


where  is  the  theoretical  linear  expression  based  on  Rytcv's  theory. 
For  conservatism  in  safety  considerations,  it  may  be  better  to  scale 
this  up  by  a  factor  of  1.28  so  that  its  maximum  is  1.6  instead  of  1.25: 


o 

m 


1.28ot 

1+0. 16o^ 


(82) 


Considering  this  maximum  in  the  standard  deviation  of  the  log  of  the 
intensity  (i.e. ,  j*1.6),  and  for  further  conservatism,  taking  a  ">o 

deviation  from  the  average  intensity,  implies: 


max 


log-1^log 
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Figure  45.  Measurement*  of  log-intensity  standard  deviation  0  versus  theoretical  (Rytov) 

Curves  indicate  various  theories  (except  in  (c)  and  points  indicate  experimental 


That  is,  the  maximum  effect  that  turbulence  can  have  (on  a  3 0  basis)  is 
to  produce  occasional  irradiance  values  that  are  approximately  five 
times  the  average  value.  This  intensity  would  be  reduced  by  molecular 
and  aerosol  extinction,  as  well  as  by  beam  divergence  effects.  The 
divergence  due  to  turbulence  appears  to  be  much  smaller  than  a  micro¬ 
radian  (for  an  example,  see  Reference  117),  and  thus  can  be  neglected 
relative  to  the  normal  optical  divergence  (which  is  not  an  atmospheric 
effect  and  was  not  considered  in  this  study). 

CONCLUSIONS 

The  detailed  conclusions  of  the  evaluations  of  the  various  theo¬ 
ries,  experiments,  and  models  discussed  in  the  previous  section  will 
be  presented  in  this  section.  The  conclusions  wilL  be  presented  in  the 
order  in  which  the  topics  were  discussed  in  the  previous  section  in 
order  to  facilitate  reference  to  the  detailed  discussion  underlying 
each.  Conclusions  of  major  import  to  the  objectives  of  the  present 
program  are  indicated  with  an  asterisk  (*). 

Atmospheric  Molecular  Concentration 

(1) *  Water  vapor  concentration  is  highly  variable  with  time  of 
day,  time  of  year,  altitude,  and  latitude,  and  can  deviate  strongly  from 
nominal  levels  provided  in  the  standard  atmospheric  models  used  in  com¬ 
puter  code  LASER. 

(2)  Some  measurements  indicate  a  smaller  water  vapor  concentra¬ 
tion  than  standard  atmospheric  models  at  altitudes  up  to  at  least  3  km. 
Thus,  the  standard  models  wou’d  not  be  safety  conservative. 

Molecular  Line  Absorption 

(1)  Molecular  line  absorption  of  laser  radiation  is  well  repre¬ 
sented,  depending  on  the  ambient  pressure,  by  one  of  the  simple  Lorentz, 
Voigt,  or  Doppler  theories,  with  the  exception  of  the  CO  and  HF  laser 
radiation. 
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(2) *  Large  inaccuracies  can  be  incurred  when  improper  frequencies 
are  used  in  the  line-by-line  calculation. 

(3)  In  the  CO  and  HF  laser  regions,  the  cutoff  frequency  should 
be  200  cm  ^ ,  and  the  exponent  on  the  (v-Vq)  term  should  be  changed  from 
2  to  1.9. 


Molecular  Continuum  Absorption 

(1) *  There  is  disagreement  between  various  data  on  water  vapor 
continuum  absorption  in  the  2400-2800  cm  ^  region. 

(2)  There  is  uncertainty  in  the  form  of  the  temperature  depen¬ 
dence  of  the  water  continuum  in  the  400-1400  cm  ^  region. 

Molecular  Scattering 

Molecular  scattering  is  well  explained  by  the  classical  Rayleigh 
express  ion. 


Aerosol  (extinction  Theory 

The  theory  for  aerosol  extinction  is  well  established  and  verified 
by  experiment. 

Aerosol  (extinction  Models 

A  number  of  models  are  available  that  use  easily  measured  param¬ 
eters  such  as  visibility,  relative  humidity,  and  wind  speed  as  inputs; 
however,  none  of  these  models  have  been  verified  to  the  extent  that  re¬ 
liable  predictions  are  possible. 

Measurement  Results 

(1) *  Extinction  due  to  aerosols  can  change  by  more  than  three 
orders  of  magnitude  during  fog  development  and  dispersal. 

(2) *  Extensive  OPAQUE  measurements  (and  others)  indicate  that 
extinction  in  the  infrared  cannot  be  reliably  predicted  from  visible 
extinction  information. 
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Turbulence  Conditions 

(1)  Low-altitude  turbulence  conditions  appear  adequately  pre¬ 
dicted  for  present  purposes  by  Hufnagel's  1978  model  [85], 

(2)  High-altitude  turbulence  conditions  are  more  difficult  to 
treat,  but  Hufnagel's  1974  model  [101]  should  be  adequate. 

Turbulence  Effects 

(1)  The  linear  effects  of  turbulerce  on  the  intensity  variance 
is  well  represented  by  Rytov's  theory  [106], 

(2) *  The  turbulence  effects  become  saturated,  so  the  maximum  ex¬ 
pected  intensity  (on  a  3 o  basis)  is  less  than  a  factor  of  five  larger 
than  the  average  intensity. 

RECOMMENDATIONS 

Recommendations  in  this  section  are  based  on  the  evaluations  and 
conclusions  discussed  in  previous  sections.  Although  the  recommenda¬ 
tions  will  be  treated  in  the  order  in  which  specific  topics  were  dis¬ 
cussed  earlier,  they  do  not  directly  correspond,  on  a  one-to-one  basis, 
to  the  conclusions.  Since  the  bases  for  these  recommendations  are  dis¬ 
cussed  in  detail  In  the  "Evaluation"  section,  the  following  listing  is 
intended  to  be  concise. 

Atmospheric  Molecular  Concentration 

Whenever  possible,  the  safety  standards  should  be  based  on  current 
local  atmospheric  concentrations  (especially  with  regard  to  humidity) 
rather  than  relying  on  standard  atmospheric  models. 

Molecular  Line  Absorption 

(1)  The  AFGL  LASER  code  (with  inputs  described  below)  should  be 
used  to  evaluate  the  molecular  line  absorption  at  several  parametric 
conditions  of  pressure,  temperature,  and  humidity  for  the  laser  lines  of 
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Interest,  to  be  used  in  semiempirical  equations  described  in  the  next 
item.  Inputs  to  the  code  should  be: 

•  The  AFCL  Line  Parameter  Compilation  updated  with  the 
latest  experimental  information. 

•  The  best  available  information  on  laser  frequencies. 

•  The  best  available  information  on  the  proper  line  shape 
and  cutoff  frequency  (i.e.,  the  maximum  distance  from 
laser  line  that  absorption  lines  are  considered)  for  the 
frequency  range  of  interest  (especially  in  the  CO  and  HF 
laser  regions). 

(2)  A  simple  algorithm  based  on  the  parametric  LASER  code  re¬ 
sults  should  be  used  to  establish  safety  standards  for  the  actual  atmo¬ 
spheric  path  of  concern. 

Molecular  Continuum  Absorption 

(1)  The  molecular  continuum  models  incorporated  in  the  AFGL  LASER 
code  (i.e.,  H^O,  N2,  C02)  should  he  used  in  the  evaluation  recommended 
under  Molecular  Line  Absorption. 

(2)  For  conservatism  in  safety  considerations,  all  other  continua 
should  be  neglected. 

(3)  Efforts  should  be  made  to  incorporate  Kunde  and  Maguire’s 
temperature  functions  [62]  for  the  400-1400-cm  ^  water  continuum  into 
the  AFGL  model. 

Molecular  Scattering 

The  existing  molecular  scattering  model  incorporated  in  the  AFGL 
LASER  cede  should  be  used  in  the  evaluation  recommended  under  Molecular 
Line  Absorption. 

Aerosol  Extinction 

Because  of  large  uncertainties  in  the  results  of  current  models, 
the  safety-conservative  Clear  Model  of  LASER  should  be  used,  unless  de¬ 
tailed  aerosol  size  distributions  and  composition  are  available. 
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Turbulence  Conditions 


(1)  Hufnagel's  1978  [85]  low-altitude  atmospheric  turbulence 
model  for  maximum  insolation  should  be  used  at  altitudes  below  the  first 
strong  inversion. 

(2)  Hufnagel's  1974  [101]  high- altitude  atmospheric  turbulence 
model,  with  a  safety  factor  of  1.1,  should  be  used  at  altitudes  above 
the  first  strong  inversion. 


Turbulence  Effects 

The  empirical  model  of  Johnson,  et  al  f 109 J  (scaled  up  by  a  fac¬ 
tor  of  1.2,  Eq.  82)  should  be  used  with  Rytov's  linear  theory  [107,  and 
Eq.  79J . 

Additional  Recommendations 

The  actual  form  of  the  model  that  should  be  used  in  setting  safe¬ 
ty  stand. rds  will  now  be  discussed  briefly.  It  is  assumed  that  this 
overall  model  should  be  highly  user  oriented  so  chat  people  not  familiar 
with  the  details  of  the  effects  can  easily  and  rapidly  obtain  useful 
results.  Due  to  the  complex  nature  of  the  problem,  a  computer  code  is 
recommended.  Ancillary  charts  and  monographs  should  also  be  prepared 
to  estimate  some  of  the  simpler  parts  of  the  problem  (e.g. ,  turbulence 
effects  and  molecular  continuum  absorption).  On  the  other  hand,  the 
code  should  not  be  required  to  evaluate  line-by-line  absorptl  ecause 
of  the  complexity  and  extensive  computer  time  involved.  These  w„lcula- 
tions  would  be  parameterized  and  fit  with  a  simple  analytic  function 
involving  temperature,  pressure,  and  humidity.  The  coefficients  of  this 
function,  for  the  entire  parameterized  range  of  conditions  required, 
should  be  stored  in  the  code  recommended  here  (i.e.,  simpler  than  line- 
by-line  codes).  To  treat  nonuniform,  slant,  or  vertical  paths  in  the 
atmosphere,  a  simple  but  physically  correct  interpolation  scheme,  which 
operates  on  the  absorption  coefficient,  should  be  devised. 

This  function  could  also  incorporate  the  molecular  continuum  ab¬ 
sorption  as  well  as  molecular  scattering  and  aerosol  extinction,  but 


evaluating  these  separately  would  be  more  desirable  since  they  involve 
such  simple  expressions  and,  in  some  cases,  their  individual  contribu¬ 
tions  may  be  desired.  Thus,  the  code  would  also  contain  simple  analytic 
expressions  for  molecular  continuum  absorption  and  molecular  scattering. 
Aerosol  extinction  requires  a  table  of  spectral  coefficients  in  addition 
to  a  simple  analytic  expression,  but  the  tables  are  not  large.  Finally, 
analytic  expressions  would  be  incorporated  to  estimate  the  effects  of 
atmospheric  turbulence.  These  should  employ  Hufnagel's  low-  or  high- 
altitude  model  [85,101]  (or  both,  depending  on  the  path)  to  estimate  the 
level  of  atmospheric  turbulence,  and  use  Rytov's  linear  expression  [107] 
or  Johnson's  empirical  expression  [109]  for  saturated  conditions.  Non- 
uniform  paths  should  be  subdivided  into  a  specified  number  of  path  ele¬ 
ments  that  would  then  be  considered  uniform. 

Several  different  options  should  be  made  available  to  the  user  so 
that  information  with  as  much  detail  as  is  available  can  be  used  in  the 
calculations.  For  exampie,  if  actual  measurements  or  other  available  „ 
information  indicates  atmospheric  molecular  concentrations  different 
than  nominal  values,  then  these  better  talues  could  be  input.  Similar 
alternate  options  should  be  available  for  measured  aerosol  size  distri¬ 
butions,  but  such  calculations  should  be  limited  to  approximate  expres¬ 
sions  for  the  extinction  efficiency  factor  rather  than  the  complete  Mie 
calculations.  Some  lasers  emit  at  several  different  frequencies,  and 
the  power  distribution  between  the  various  transitions  should  be  speci¬ 
fiable.  If  such  complete  information  is  not  available,  nominal  distri¬ 
butions  should  be  stored  within  the  code.  The  output  of  the  code  should 
also  be  straightforward  and  fully  labeled  so  that  it  is  easy  to  under¬ 
stand. 

It  is  further  recommended  that  this  code  have  the  capability  of 
evaluating  the  probability  that  safe  irradiance  levels  would  not  be  ex¬ 
ceeded,  at  specified  ranges,  as  a  function  cf  tune  of  year,  by  using 
statistics  on  meteorology,  atmospheric  constituents  (molecular  and  aero¬ 
sol),  and  turbulence.  This  feature  should  be  available  in  operational 
and  conceptual  test  planning. 
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APPENDIX  A.  MOLECULAR  CONTINUUM  ABSORPTION  KEAS 
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with  increasing  temperatura  much 
faster  than  theory  predicts. 
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